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ABSTRACT 


Propagation of long waves in a two dimensional systems is 
analysed. An acoustic approximation is used and "boundary con- 
ditions simulating typical industrial piping systems are formu- 
lated. The thesis consists of three main parts. In the first 
part a general solution of the boundary value problem is obtained. 
In the second part series solutions of the equations of motion are 
generated, and a physical interpretation of results is given. A 
numerical application and a discussion of the problem are given in 
the last part. 
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INTRODUCTION 


The characteristics of motion of a fluid, whether steady or 
unsteady, depend on the nature of the fluid, the imposed initial 
conditions and on the nature of the system in which the motion takes 
place. In the large class of unsteady motions ,of compressible 
fluids, the periodic motion or wave propagation is particularly im- 
portant from the point of view of both pure and applied sciences and 
has been the subject of extensive analysis. However, except for 
essentially sircple systems the motion of pressiire waves is very com- 
plex. Generally, the waves may be reflected, diffracted and dis- 
torted. They may be subject to attenuation and decay or to inhomo- 
geneous distribution. The more the symmetry of a system departs 
from a simple rectangular, cylindrical or spherical geometry and the 
more the initial conditions differ from that appropriate to the 
boundary, the more these effects will be pronounced. 

Besides a purely scientific interest in the problem of propa- 
gation of waves, there are many engineering reasons for a better 
understanding of wave motion. Of theoretical and practical interest 
is the determination of the basic modes of motion and their phase 
and group velocities. Also knowledge of the resonant and the cut-off 
freq.uencies as well as the determination of the reflection and trans- 
mission coefficients are of substantial interest. 

f 

Propagation of waves in curved ducts and pipes belong to the 
class of motion which is characterized by wave patterns totally 
different from those known in straight ducts or in iinlimited space. 


1 
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The curvilinear. boundaries are -responsible for the appearance of a 
continuous steinding radial wave which in turn affects the trans- 
mitted tangential waves . The propagation in curvilinear waveguides 
or bent ducts is difficult to anal/ze, and mathematical models de- 
veloped to date are complicated. Perhaps that is one of the reasons 
why relatively few papers are available on this subject. 

The purpose of this treatese is to solve the problem of propa- 
gation of waves in curved ducts. It is xntended to obtain a solu- 
tion for long acoustic waves in sli^tly and sharply bend ducts. 

This problem has been only partially analyzed by varioxis authors . 

The first recognition of the problem of propagation of pressure 
=1 

waves in a curved conduit, as distinct from the motion in a straight 
line, was formulated by Rayleigh In a short, brilliant expose 
he demonstrated that long waves -in. a curved pipe of infinitesimal 
cross section behave exactly as in a straight pipe. The curvature 
of the pipe is of no importance. His analysis is based on the lin- 
earized equation of motion and conclusions are not valid if the 
cross section of the pipe is finite. However, what is important is 
that he establishes a limit to which equations describing wave motion 
in pipes of arbitrary shape must tend. 

Subsequent research of wave motion in curvilinear ducts is al- 
most exclusively anal;^ieall . The majority of papers treat the prop- 
agation of electromagnetic waves in curved ducts. Only a few discuss 
the propagation of sound waves. Interestingly enough, along with 
analytical fonaulation of the behavior of waves in bends there 
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appeared a series of papers dealing with mathematics needed to solve 
the physical problem of waveguides . This parallel effort indicates 
that the solution of the problem of bent waveguides requires mathe- 
matical formulations and techniques not generally ^available . 

To review the literature on the wave propagation in curved 
ducts it is desirable to discuss the work done by grouping studies 
by adopted boundary conditions. In several instances studies orig- 
inally intended to be general were actually reduced to some specific, 
restricted cases by simplification of boundary conditions -in order 

to solve the problem. 

(2 1 ( 3 ) 

In 1905 Kalahne' ' and in 19^6 Jouguet' and a decade later 

(J+) 

Karpman considered oscillations in rigid toroids . Jouguet 
approaches the problem by the method of perturbations while Kalahne 
and Karpman use the method of separation of variables. Both Jouguet 
and Karpman intended to study propagation but simplified the problem 
by adopting closed systems. A similar standing wave problem, using 
a configuration of an elbow connected to two segments of straight 
tube, was treated in 1930 by Matchlnskl^^^ . A very extensive and 
complete treatment of the behavior of electromagnetic standing waves 
in circular cross sections of simple or coaxial tubes is given in 
the text ''Sleet romagnetic Waveguides and Cavities" by Goubau^^^ 

As interesting as their developments are they contribute little 
to the determination of parameters of progressing waves. However, 
a number of papers published between 1939 and 1969 did .contribute 
to the problem., Buchholz*^^^ , using .separation of variables 
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obtains a solution for propagation of electromagnetic waves in 
slightly bent waveguides of infinite length. He calculates an ex- 
pression for a wave propagation constant and draws several general 
conclusions about behavior 'of waves in bends. Buchholz’s paper is 
the first contribution to the problem of progressive waves. The 

problem of the infinitely long bend was also treated by Krasnushkin^ \ 

(9) flO^ 

Grigor’yan , and by Victor ov and Zubova^ \ Krasnushkin ap- 
proaches the problem by the method of separation of variables but 
in view of mathematical difficulties proposes a perturbation method 
and treats the simplified case of the slightly bent tubes. Grigor'yai^ 
solves the differential equation by expansion of the cross-product of 
Bessel functions into a Taylor series. He tries to obtain an algor- 
ithm of sufficient generality to be applicable to all possible 
impedances of the waveguide walls. Coefficients of his series con- 
sists of a sequence of Wronskians as determined by Basset The 
method is only partially successful. Grigor’yan obtains correct 
general information on amplitude and distribution of the radial os- 
cillations but his basic mode wave number does not verify the dif- 
ferential equation except for the Raleigh case of very narrow pipe. 
Victorov and Zubova treat short wave propagation in a solid 

layer and outline a possible solution for an infinite coil of larg,e 

( 12 ) 

radius of curvature. Finally, Voskresenskin and Sveshnikov et 

(13) 

al. approach the problem of infinite bend using the separation 

of variables. However, like several other authors, they do not 
solve their equations but limit themselves to general discussion. 
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The contrihutions to the solution of wave propagation in curved 

ducts which were reviewed considered only one aspect of. the problem, 

the functions which depend on radii and on boundary conditions of 

the bend walls . The dependence on the longitudinal coordinate , that 

is, the angular propagation was greatly simplified by assuming an 

endless coil. However, two papers consider the angular dependence 

(l4) 

in detail. Rice , using matrix equations establishes the basic 
mode propagation constant for electromagnetic waves and calculates 
the coefficient of reflection of waves incident upon a bend. Bas- 
ically, he treats a junction between a straight and a curved cylin- 
drical duct. His procedure - the theory of matrices - is an obvious 
choice for a system capable of supporting several modes of wave 
motion. Arlinger^^^^ , using the method of perturbation, presents 
a problem of compressible potential flow in a bend of a duet of cir- 
cular cross section. Although this is not a case of propagation of 
waves, Arlinger's general solution for potential motion in a bend 
and at a junction is closely related to our problem. The common 
feat\ire of the two papers is need for extensive supporting mathe- 
matics. Both authors obtain solutions which contain definite in- 
tegrals which cannot be integrated by ordinary means . To by-pass 
this technical difficulty Rice supplies tabulated values in an 
appendix. Arlinger resorts to integration by Simpson's rule. 

A number of papers on roots of equations containing cylindrical 


functions.,' (cross-products of Bessel functions) were found in the 
process of verifying the existing literature on the wave motion in 
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cylindrical coordinates. Unfortunately most of these papers, (by 
McMahon^j*"^^ Truell^^'^^ Dwi^t^^^^ and Kline^^^^) are limited to 
research of characteristic values of the argument of the Bessel 
functions. Therefore, their usefulness is limited to problems con- 
cerned with finding resonant conditions in closed systems . However, 
one monograph by Buckens^^*^^ , does apply directly to -relationships 
for evaluating the order of cross paroducts of Bessel- functions of 
real and of imaginary order, and was used in developing the ex- 
pressions for decaying oscillations. 

Among the few treatises on propagation of sound in bends are 

(21 22) 

those by Miles ’ in -sdiich he establishes an analogy between 

propagation of sound and an electrical transmission line are most 

inportant. The method is then applied to a right angle joint of 

rectangular tubes. Actually, Miles is considering a mitered 90 ° 

bend with the outside corner chamfered. "This greatly simplifies his 

( 5 ) 

analysis . Bends of this type are discussed by Goubau and of 
course constitute a special case of "reflectionless" bends which may 
be used for a selective range of frequencies. 

(23-25) 

The work of Miles was successfully checked by Lippert 
Lippert presents an experimental study of. sound wave propagation in 
mitered bends of various angles . The experiments were conducted o ■ 
over a wide range .of frequencies and show that long waves in mitered 
bends propagate with insignificant reflections . In- an interesting 
experiment with. a rounded 90°»bend (the outside wall of the miter 
joint of the duct was made circular with a radius equal to the duct 
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width) Lippert shows that the transmissivity was far superior to that 
of an ordinary miter joint over a very wide range of frequencies. 
Finally,' it should be mentioned that Morse and Feshbach^^^^ propose 
an approximate method for the solution of the problem of long wave 
propagation in a two-dimensional 90° miter -bend. A series solution 
of the wave equation is written for the two infinite branches of the 
two-dimensional . duct . The square area of the elbow proper is anal- 
yzed by the method of conformal transformation. 

An analytical and experimental program on wave propagation in 

(27) 

flexible pipes was carried out by HASA . Propagation of small 
sinusoidal perturbations was studied in a line including a 90 ° elbow 
with' flexible supports. The experiments were supported by an analy- 
sis based on the acoustic impedance theory which described the dy- 
namic behavior of the system. The conclusion of the experiments is 
that the elbow created no noticeable effects. 

Propagation of waves in bends and more particularly in circular 
bends still requires fundamental analysis. The, thus far, obtained 
algorithms for computation of the wave motion are very cumbersome 
and the past research covers only certain aspects of the wave motion 
in curved bends. 

The general object of this dissertation is to analyze the be- 
havior of long waves in curved conduits. The purpose of the analysis 
is to determine the distribution and variation of the vibrational 
velocities, and the phase of motion, at any point in the systems 
under consideration. The particular problem is to determine a 
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mathematical model of the 'wave motion that will yield, ■vd.th reason- 
able accuracy, meaningful answers and allow a physical interpreta- 
tion of results 

To satisfy these reg,uirements the propagation of acoustic waves 
was selected for model. The wave transmission through bends of ar- 
bitrary sharpness is solved, for the first time, as boundary value 
problem using the characteristic functions. 

f ] 

The following conditions are assumed to exist. There is a sus- 
tained, continuous and steady propagation of long acoustic waves in 
a two-dimensional curved channel. Long acoustic waves are under- 
stood to be harmonic vibra-tions of infinitely small amplitude of a 
compressible, in-viscid fl-uid. The wavelength is at least two orders 
of magnitude larger than the width of the channel. The pressure 
distribution -sriLthin the channel satisfies the wave eq^uation. The 
radius of curvature of the bends is unrestricted. 

In this study two acoustic systems are considered which allow 
determination of the basic modes of motion and describe the transi- 
tion and distortion of plane wa'ves as they propagate do'wn the curved 
channel. A detailed study of the beha-vLor of waves in junctions be- 
tween strai^t and curved ducts iTill be given. 

The mathematical treatment of the problem utilizes the method 
of separation of variables.. Solutions and expressions for principal 
modes of the wave are obtained by using the linearized equation of 
motion solved for its characteristic values. 

This original approach required a novel use of Bessel functions 
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to determine the characteristic values of the steady and the decay- 
ing fields of motion. 

The unavoidable approximations in the numerical solutions of 
the cylindrical functions are due to use of series expansion of Bes- 
sel functions and from restrictions necessary to solve infinite ma- 
trices. 
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1.0 GEUTEML SOLUTIOI 
1.1 Formulation of the Problem 

1.10 Equation of motion . - Bee axis e the complexity of the problem 
considered requires a simplified mathematical model-, a linear equa- 
tion of motion, the wave equation. 




( 1 . 1 ) 


will be used. It is generally tractable and should yield resxilts 
approximating, to a degree, the real phenomena in curved pipes. 
Equation (l.l) is well known and there is no necessity to derive it. 
It is, however, important to indicate here the usually omitted dis- 
cussion of the degree of approximation resulting from use of a lin- 
earized equation and the nature of applied assixmptions . The best 
way to present all approximations and assumptions is to derive the 
wave equation from a most general equation of motion. 

I 

The derivation given hereafter is classical. Elements of it may 

(28) ( 29 ) 

be found in most'Of ithe advanced texts , . The most complete 

■ (2(i) 

derivation is given- by Morse and Feshbach' . 

The sxxm of forces on a fluid element equals the acceleration 
of the element, DV/Dt , times its mass p dx dy dz. 

The very general equations of motion, without restrictions, as 
to the compressibility, viscous forces, and turbulence, is 

P ^=-- V(P + ¥) + V .7 (1.2) 

DV 

err is the substantial derivative of velocity vector, a sum of the 
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local contribution to acceleration and the convective con- 

o t 

tritution due to translation 

^ = V ■ ~ V{V • V) - Y X curl V = ^ W - V x curl V 

The -VW = F- is the external potential of conservative forces of 

gravity and V • x is the viscous force on element per unit volume 

/ 

but because of the tensorial nature of x, V * x is not a simple 
divergence. 

To simplify the problem we restrict our analysis to small aaip- 
litude displacements and velocities . ¥e also assvime that the com- 
pressible fluid under consideration exhibits no free surface, is 

homogeneous throughout the domain. In this case all terms in equa- 

2 

tion ( 1 . 2 ) involving V can be neglected and^ assigning that the 
compressional viscosity -will be small, the influence of viscosity 
will be negligible. Furthermore, in the equation of motion of a 
homogeneous fluid, in the absence of a free surface, the forces due 
to the external potential (gravity) can be disregarded if pressure 
P is the excess pressure p over that of the fluid at rest, P^. 
With these assxunptions , the equation of motion reduces to 


P 



( 1 . 3 ) 


IFow we may write V grad ^ where 4i(x,y,a,t) is a scalar function 
or potential whose partial derivatives with respect to x,y,z are 
the components of Velocity Y in those directions. Substituting, 
differentiating equation (1.3) with respect to time and since 
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= - — div V as established in appendix 1, equation (Al.l), 
dx X 

the eq,uation of motion for small oscillations, of a compressible, 
homogeneous and inviseid fluid at rest is 


2a - 1 


V^{j) = — 2" 


at-^ 


(l.l) 


where = l/xp* This is the linearized wave equation in tei*ms of 
the velocity potential. It will he convenient to express equation 
(l.l) in polar coordinates to fit the geometry of circular curved 
pipes . 

The two-dimensional equation of motion of the analytical model 
is then 


£ 1+ 1 ^ (1 4 ) 

3r2 r^ 302 c^ at^ 

The particle velocities (the vibrational velocities) are given by 

v(r,0,t) = “ (the tangential component) 
r do 

u(r,e,t) = (the radial component) 

and the excess pressure by 
p(r,0,t) == “P H" _ 

1.T1 Boundaiiy conditions and physical systems to be considered . - The 
boundary conditions which the solution must satisfy depend on the 
physical configuration of the system to be analyzed. On the othr 
hand we are not entirely free to choose any boundary conditions 
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because the solution of the hyperbolic equation to have a physical 
meaning must be done over an open boundary at least in ’one space 
dimension which may be along the time axis . In the case of the prop- 
agation of waves in pipes the walls are a closed boundary, the origin 
and the far end sections may be opened or closed. The time axis is 
usually open. On the walls and on the end sections the Diriehlet 
and the Neumann conditions may be used. The Diriehlet condition re- 
quires that the potential be specified on the boundary. The Neumann 
condition requires that the normal component of the velocity on the 
boundary be specified. 

Two physical systems shown on figure 1 will be considered. 

Both approximate typical industrial piping. The first consists of 
a rigid infinite bend approximating a coil. At the inlet section 
of the coil pulsations are generated by a hypothetical piston of 
infinitesimal displacements . Conditions at the far end section of 
nhe bend are of no consequence because assumption of an infinite 
coil implicitly states that the far end section condition will not 
contribute to the solution; only a simple wave in a positive 6 di- 
rection will be considered. 

The second system consists of a circular bend followed by a 
straight, infinite duct. As in the case of the infinite coil, at 
the inlet section of the bend pulsations are generated by a hypo- 
thetical piston of infinitesimal displacement. The far end of the 
system will be at infinity. The 'system of coordinates and nomen- 
clature used are shown on figure 1. 







15 


It will be assumed that the walls of the curved and the strai^t 

ducts will he perfectly rigid so that the Newman condition will 

apply. Generally, walls may have an admittance (reciprocal of the 

normal impedance) different from zero and a very general solution 

may he obtained on such a basis. An extensive discussion of this 

aspect of the boundary conditions and formulation of this problem 

(9) 

may be found in Grigor’yan's paper . Our assumption of perfectly 
stiff walls and hence, the requirement that, at the walls, the nor- 
mal component of the vibrational velocity vanishes is a very good 
approximation of actual acoustic systems . 

1,2. Determination of the Characteristic Values and of the 
Characteristic Functions 

The linearized wave equation in cylindrical coordinates is 
known to be separable in coordinates proper for the boundary. Equa- 
tion (1.4) may be broken up into a set of ordinary differential 
equations, each including a separation constant. Consequently, no 
other method of solving our boundary value problem will be contem- 
plated in spite of the fact that several approximate methods, such 
as the method of perturbation, have been found useful and leading 
to meaningful, although generally, restricted results. 

To solve equation (l.U), we assume a solution of -the form 
= R(r)e(0)T(t) 

By separation of variables we have 

IT" ,2 j m i(kct+a) i(mt+a) 

^ = -k^ and T = e " = e 
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This solution assumes that the dependence on time is simple harmonic 

(•with a possible phase lag a) and that equation is, in fact, 

t 

a scalar Helmholtz equation. 

0 " 

Next, Q — = The solution of the angle depending term may 

I 

be periodic or linear, depending on the value of the separation con- 
stant . 

0 = a cos v0 + b sin v6 v 0 

V V 

9 = c0 + d \)s=0 

Finally, E" + |- R' + ^ E = 0 

and 

R = A J (kr) + B Y (kr) 

V V V V 

■which is the characteristic function of the problem. 

Superposition of solutions is allowed^ for -equation (l.U) is 
linear and a general solution may be written in the form 




i(mt+a) 


(a cos v0 sin v0)rA J (kr) + B Y (kr)"]dv 
V L\)v vv -J. 


+ + d)[AQjQ(kr) + BQYQ>(kr)] 


(1.5) 


where C is a set of points in the complex plane, to be determined 
in order to satis:^ the boundary conditions . 

The solution for v = 0 might have been included in the in- 
tegral terms but in order to show more explicitly the linear depend- 
ence on 0 of this solution, it has been written separately. 
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Before proceeding with application of the hoimdary conditions 
we will consider the well known technique of splitting the potential 
field into the incident and the scattered waves 

(fc = A + (f> 

The incident plane waves = v^ exp i(it)t - kr sin 0) propagate 
from the vibrating pistons into the bend until they strike the curved 
wall facing the piston. However, shortly after leaving the piston, 
the waves cease to be plane because they mil be strongly diffracted 
on the curved inner wall of the bend. It is noted here that the 
work on diffraction of long waves by the edge of a semi-finite screen 
by Lamb^^^^ in 1906, and generalized by Penny and Price in 1952^^*^^, 
is not valid on ciirved walls. 

The reflected, scattered waves which will result from the im- 
pact on the outer curved wall, will form cusped fronts which have 

(31) 

been experimentally studied by Wood , These fronts will carzy 
the energy of original waves down the bend. As the front progresses 
down the bend it will spread, and assume some new form. The initial 
condition of the reflected wave is that of the diffracted wave 
reaching the outer curved wall. Because the degree of diffraction 
of the plane waves is unknown, the initial conditions of the re- 
flected wave are undeterminable . In view of the difficulty of ob- 
taining proper expressions for the incident and the reflected waves, 
this approach is not applicable. Consequently, the case of motion 
in bent ducts (which lack axial symmetry leading to the diffraction 
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of waves and multiple, endless reflections) requires a solution which 
would not include the identification of the component waves. 

To satisfy the partial differential equation and the boundary- 
conditions for perfectly rigid circular walls, a characteristic 
equation will be foimd whose roots will be the characteristic values 

t 

of the problem: a set of values of the separation parameter v 

which will yield a nontrivial solution of the problem. Differenti- 
ating equation (1.3) with respect to r and equating to zero we 
obtain for the two circular boundaries 

A J'(kR, ) + B Y'(kRj = 0 

1 V # 0 (1.6a) 




(1.6b) 


Conditions equation (l.6b) give Bq = -A 
and J^(kR^)Y^^(kR 2 ) - J^{kR 2 )Yj^(kE^) = 0 


J^CkR^^) 

0 Y^(kE^) 


Since jQ(kr) = J^(kr) and YQ’(kr) = - Y^(kr), equation 
be written 


(1.7a) 


(L.'7a) can 


J^(kR^)Y^(kR2) - J^(kR2)Y^(kR^) = 0 


(1.7b) 


For the case of (kr) greater than the first root of Y^(kR^) the 
cross-product (eq. (1.7b)) will yield a series of solutions. In the 
present case where 0 < (kr) < (kEg) Ij with the steep slope of 
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Y^(kr) and relatively moderate increase of J^(kr) over the range 
(kR^) to (kEg), there is no solution for situation (1.7b) in the 
range under consideration. 

Consequently, the solution v = 0 cannot be considered. Equa- 

t 

tions (1.6a) give 

~ \ Y'(kR, ) 

V 1 

and J^(kEj_)Y^.(kRg) - J^(kR 2 )Y^(kE^) = 0 (1.8) 

Considering R^ = aR^^ ("a" is a proportionality constant) and 
using the relation 

cos(ttv) J'(kr) - J’(kr) 

Y^(kr) - — — 

sin('irv ) 

equation (1.8) simplifies to 

There are no known tables for the roots v of the character- 

m 

istic equation (1.9). To evaluate the 'W'® expend the and 

in terms of increasing powers of the argument (akR^) and 
(kE^). If we limit the expansion in the first approximation to the 
first term: 

J*(kr) (kr)'^ ^ + . . . 

^ 2^'r(v + 1) 

J' (kR) - (kr)“^"^ 

2"'’ (1 - v) 


+ . . 
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and substitute into equation (2.5) we obtain 




= 0 


and finally 


V ( v-1 -v-1^ 


( 1 . 10 ) 


Solution V = 0 has been already rejected. Therefore, the 
only acceptable solution must satisfy the equation 


v-1 -v-1 

a =s a 


which may be put in the form 


2v 2vlna ^ 

a - 1 or e =1 


Hence, 2vlna = 2miTi 

,, . . miTx 

that IS V = 

m Ina 


(m = 1, 2, 3 


(i.ir) 


Better approximations will be given by the second and following 


terms of expansion of-^ J' and J' 

V -V 


j« ( 2 ) v(z)^~^ _ (2 + v)z^'*‘^ ^ (Ij. + v)z^'^^- 


2^r(v + 1) 2^ (2 + v) 2 ^r(3 + v)r(3) 


(i. 12 ) 




-V z 


-v-1 


(2 - v)z^~'^ , (4 - v)z^ 

“ — — _ -J. 


r(i - v)2 ^ 2^"^r(2 - v) 2^“''r(3)r(3 - v) 


Since z = kr is small for k small , we neglect here the 
third term of these expansions. 
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Equation (1.9) gives : 



(2 + v)(kR^)^‘^^ 

^ i 

-v(kaRj^) ^ (2 - v)(kaR^)^~^ 

_2^r(v + i) 

2 ^ • Ur (2 + v)_ 

_ (l - v)2"''' 2“^ ■ k *r(2 - v)_ 


(kaR^)^“^ (2 + v)(kaRj^)^‘*''' 

-v(kR^)“'^~^ 

2"^r(v + 1) 2^ • Ur(2 + v) „ 

_r(i - v)2“'^ 


(2 


it • 

2"^r(2 - v)_ 


where , 

r(v + i) = vr(v) 


r(i - v) 


IT 1 

sih trv r(v) 


(1.13)' 


r(2 + v) = (i + v)r(i + v) = v(i + v)r(v) 


r(2 - v) - (1 - v)r(l - v) = (1 - v)ir/sin(7rv) • r(v) 
r(i - v) = -vr(-v) 

r(l + m + v) = (l + v)(2 + v) . . . (m + v)r(l + v) 

After algebraic manipulations and dropping the term containing 
2 

the (kRj^) we have 


sin irv 


l6irv(l + v)’(l - v) 


1+ v(l - v )(2 + v)^a ^ ^ - 


Hv2(l - v)(l + v)(kE^)“^(a^"^ 


a^'^^]+ (1 + v)(2 - v) 

- = 0 (l.lA) 
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In order to improve on the first approximation, (eq^. (1.11)), 
■we now assixme 


V 


miri - . 

in* 


ln(l + e ) 
m 

In a 


(1.15) 


where is a small quantity. 

Solving and remembering that 

In a = ln|zj + i(argz + 2nir) n = 0, ± 1, . . . 


we have 

' Ur, ^ -unri-lna rnTfi+lna') 

0 = (1 - v)(2 + v)^e - e J 

+ (1 + v)(2 - 
1 - l,v(l - v)(l + 

0 = (1 - v)(2 + v)(-l)“(a-^ + a'^^) + (1 + v)(2 - v)(-l)“(a"^ - a"^^) 

•. l,v(l - v)(KR,)-"(-l)-‘ 

since — ; — = 1 - e 

1 + e 

(2 - v^)(a"^ - a"^^) + 2v(l - v2)(kR^)“^ ° 

(2 - v2)(kR^)^ 

= “ 4v(l - v'^) _ 


E 


m 


“ i 



(a2 - 1) 

2 + 

2 2' 
m^TT^ ) 

(in a)^_ 

4mTT 

~ 1 n r> 

m^TT^ 

» f 

^ Tin 


In a 


(l.l6 ) 
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and 


_ . J mfr 
V = 

m tin a 




(kEj_)2(a,2 ^ 1) 


2 + 


m^TT^ 
(in. a)^ 


iiniTr 


1 + 


2 2 

(in a)' 


(1.17) 


For m = 0 


ln(l + £ ) 
o 

V = 

o In a 


m = 1, 2, 3 


Substituting into equation (1*14) and neglecting terns in we 


have 


2 - 


ln(l + 

In a / 


— r ~ a(l + £ ) 

a(l-+ £ ) o 

o 


2 + 


ln(l + £ ) 
o 


In a 


1 + £ 


a 


1 + £ 


+ ij- 


ln(l + £ ) 


O' 


In a 




-2 P ^o 


a(l + e ) 

O y 


= 0 


(1.18) 


Solving, we have 
/ ln(l + £ )A 1 - £ 


2 - 


In a 


f ln(l + £ )\l + £ 

° a(l + e). + 2 + = 2_) 2. 

V In a / a 


a(l - e) + U 


ln(l + £ ) 


2e 


T (kR, )■ - — = 0 

In a 1 a 


-K" - 1) 


+ 2 


ln(l + £^) 

In a "^o 


(-1) 


8e ln(l + £ ) „ 

+ — 2. 2_ (kjj = 0 

a In a 1 


Finally 


=o * ®o) 

in a ^ 


- 1) 


- (kK )~^ + a + - 
a- 1 a 
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ln(l + e) = e 

_ 2(a^ - l) 

In a 4(kR^)‘'^ + a^ + 1 


and 


V 

o 


2 ( 


a‘ 


l)/ln 


a 


4(kR^)“^ + a^ 


+ 1 


(1.19) 


Since the result depends on we retrace our steps and solve 

equation (1.18) without neglecting terns in v^. After algebraic 
manipulations we obtain a new expression for 


V ■ 

o 


2(a^ - l)/ln a 
4(kK ) ^ +'a^ + 1 + 


- 1 

In a 


( 1 . 20 ) 


There is an infinite set of pure imaginary .roots given by equa- 

tion (1*17) and one sin^e real root given by equation (1.20). The 
uniqueness of the obtained real root remains to be verified. Ex- 


panding equation (1-9) by means of series (1-12) and substituting for 

V v^(ln a)^ 

a the power series a^ ^ 1 + v In 'a + ^ ... we obtain, 

in the first approximation. 



= (kE^) 


4(kE^) 


a^ - 1 
2 In a 


This result was obtained with only first two terms of the series 
for a'’. When three terms of this series are used and when small 
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of the fourth order are neglected, we obtain 


- 1 

In a 


4(kK )”‘^ + + 1 + ' V — - - - (a^ 

1 In a 


- l) In 


Equation (l.20) is thus Terified by series expansion and uniqueness 


of the root established. Adding all solutions obtained with 


each of the calculated root, equation (l.S) becomes 




cos V e + sin V 0^ A J (kr) + B Y (kr) 1 (l.2l) 
' m m V \> ' V V \ ' 

m y m m m m J 


and using the previously determined value of B in terms of A 

m m 

J' (kE, ) 

V 1 

B^ = - A^ — ■ y - , and simplifying notation by adopting 
m m 1'^ 


C 

V 

m 




we obtain 


00 


( a cos V! 0 + sin v 9 
\ V m m 

m 


-Y’ (kR, )J (kr) 

V 1 V 
m m 


(k\)Y^ (kr) 
m m - 


which still further simplifies to 
C 

CO yj — 

^ \ ~ — — — (a cos V 0 + sin u 0^ J (kr)j' (kR, ) 

/. i Sin TTV V V ro- m ) V -V 1 

m \ m / m m 

m=0 *- 


-J (kr)J' (kR, ) 

-V V 1 

m m .. 


( 1.22 ) 
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The exjiression in the square bracket is the characteristic function 

of wave motion in curved ducts. For simplicity it will be denoted 

F^(r,E^) or in an abreviated form F^. Its derivative with respect 

to the argument of Bessel functions , for the inner curved wall , that 

is, at is always equal to zero. On the outside curved wall, at 

R-, F is zero for all characteristic values v . To determine 
d V m 

values of constants we differentiate equation (l.22) with re- 

m 

spect to 6 and require that, at 6=0, the tangential vibrational 


velocity of particles be equal to the velocity of the vibrating pis- 


ton 


V 

o 


imt 

e 


1 H 

r 30 


e=o 


C V 
m m 


J' (kBj 

V 1 
m 


J (Jcr)'J' ( 

-V V ’ 

m, m 


V 

o 


(1.23 ) 

where, to simplify the notation, C is written for C 

m 

To be able to write an expression for we have to use the 

orthogonality conditions satisfied by the set F^ = (kr)j' (kR^) 

-- m L m ^m 

-J_^ (kr) (kR^) , as given in appendix 2. 

m m - 

The reSTolt is 




(1.24) 
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!The integral | — is easxly reducible to a s\m of integrable 

J ^ \ 

expressions of which one vanishes . The evaluation of this integral 
is given in appendix 3. 

The result is 



rF 


- f 2 dr = 


m 


r V 


2v ^(v) 
m 


m 


(1*25) 




9F 

V 

where F^^^ = 

V 

m 

to the argument of 


and 



' indicates differentiation with respect 


The integral 



V = 
o 


constant, can be integrated using the following expansion in Bessel 

( 32 ) 

functions, as given in "Handbook of Mathematical Functions". 



2 

k 



+ J 


v+3 


+ . 



R. 


The sums of Bessel functions of increasing order ±v^ +1, 

±v^ + 3, • • . in complex numbers, converge rapxdly and the two 
first terms should give a satisfactory solution. 

Summing up, the final expression for the velocity potential is 
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OO 



D 

m 

sin(7rv ) 
m 


cos V 0 
m 


J (kr) (kR^) 

L m ~ m 



(kr) (kR^) 
m 


) 


(1.26) 


where 



( 1 . 20 ) ' 


m 


- 


mit 
In a 


(kR^)‘^ea^ - r) 


2 -+ 





HI = Ij 2, 3. . . 


( 1 . 175 ) 


In equation (1.26) only D^, proportional to a^ in equation (1.5), 
remains undetermined. 


1.3 llypical Physical Systems 

The constants must he determined hy the remaining boundary 

conditions . These conditions depend on the configuration of the sys- 
tem under consideration. 

1.30 Infinite bend . - Suppose an infinitely long circulaS* duct for ■ 
which the closest physical example is a tightly wound coil of whxch 
the pitch is negligible compared to the radius of curvature of the 
duct. The far end boundary condition for an infinite duct is that 


no reflection of waves must be considered. 
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To determine D in equation {l.lS) "we rewrite (C sin v 6 
m - mm 

+ cos v^e)j which is a convenient expression for study of standing 


m 


waves, into a form better suited for study of propagation of waves. 

0" o 

An alternate form of the solution for ~ = - is C e 
-iv^e ® 

+ D* e . The two coefficients C and D' are not known, but 

m - mm 

may be obtained in terms of and D^. Equating the two expres- 


sions 


sin(irv ) 
m 


iv0 . 

2 ^ "" -"2 


+iv 0 
m 


where D is the unknown to be determined, 
m 

The new, alternate, expression for the velocity potential is 

F 


4, = 


itot 1 
® 2 


m ==0 


fiC + D )e “ + f-iC + D 
^ m my ^ ™ 


m 


sin(7rv ) 
m 


(1.27) 


where F^ is the radius depending terra of the velocity potential, 
m 

the characteristic function of the problem. 

Let us consider the imaginary roots v - ip , with y >0, 

m m’ '^m ’ 

(m = 1, 2, 3, • . . ) . We have now: 


CO 



m=l 


(-iC -K D ) 
\ m ra/ 


-p 0 

m 


F 

V 

ra 

sin(7TV ) 
ra 


As 9 *> “, to avoid values which increase without limit with increas- 
ing values of 0 we must set 
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and 



-e 


iut 


iC 


m 


3n=l. 


m 


sin(TTV ) 
m 


F 


m 


CQ 



( 1 . 28 ) 


which represents attenuated vibration of a compressible, inviseid 

fluid. Solutions pertaining to the imaginary roots vanish shortly 

after the inlet to the duet. The solution given by the real root 

V is 
o 


<J> 


o 



+ 




F 

V 

0 

sin(7rv ) 
o 


To eliminate the waves going into the negative direction of 0 , we 


put iC = D and obtain 
o o 


iC F 
o V 

*^o sin(TTV ) ^ 


i(mt-v e) ^o^v 
o o 


i ( (ot-v e ) 
o 


sin(TTV ) 
o 


( 1 . 29 )‘ 


tdiich is a wave equation representing free and undamped oscillations 
of a congress ible , inviseid fluid. 

Adding the obtained solutions the velocity potential of infini 
tesimal waves in an infinite two-dimensional curved duct is 


(fi(r,e,t) 


D F 
o V 

o 

sin(irv ) 
o 


i ( o)t-v 0 ) 
o 


CO 



m-1 


D F 
m V 

m 

sin(iTV ) 
m 



( 1 . 30 ) 


The amplitude of the moving wave is given by the constant and 

by the characteristic function F (r,E^ ) of the first term of the 

o 
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eq.uation (1.30). The exponential function exp i(cut - v^0) describes 
the harmonic motion. The second term ’of this equation represents 
the decaying oscillations. The amplitude of these oscillations de- 
pends on the radius. Basically, the motion in curved bends is char- 
acterized by non-plane oscillations. The detailed evaluation of 

I 

these waves and oscillations will give the degree of this non- 
uniformity . 

The harmonic function exp i(o)t - depends on the charac- 

teristic value (eq. (1.20)) which is the wave number of motion 
when expressed in cylindrical coordinates. It has been named by 

/ O \ ' 

Krasnushkin^ "angular wave number"; is a nondimens ional, real 
and fractional number. 

The propagation of waves in a curved duct is profoundly influ- 
enced by the boundary conditions (i.e., by the curvature of the duct). 
Consider exp i(mt - v^e) where u, is proportional to the forcing 
frequency of the piston and the phase velocity 0 = 
ponential function may be written exp i(ajt - v^s/r) where s = r0 
is the arc length. The wave constant on every circular path is 
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It may be shown that for the case of a sli^tly bent duct the 
dereloped equations for wave propagation in bends are , in the limit , 
exactly those for a straight duct as demonstrated by jRayleigh. This 
will verify the result of the present analysis. Equation (l.20) 



gives V = A I 

° H- 1 

I -L» 


while the apparent phase velocity on a circle becomes s = y 

o 

“ — > ^ = cs velocity of sound in unbounded space. 

The angxilar wave number for short waves obtained by Buchholz' * 

(lij) ■ 

and Bice , does not meet this^condition. Using our nomenclature 
they have 

■ - * 

. / ? T . (kHj^(a - 1)^ 

Vq = a + i + I — 1 + F(a - l) 


which for a 1"* tends to - Y(kR^)^ + ^ 

(9) 

The expression for angular wave number obtained by Grigor'y^^ is 


V 

o 


a(kR^) 



It tends to .a 1, as required. However, for 
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a > 1 its value does not agree ■with, equation (l.20). In vie'w of 
the fact that Grigor'yan uses an approximate solution (expansion in 
Taylor series) while equation (l.20) is a characteristic value of 
an exact solution, a verification of the solution at the lower limit 

I 

and a progressive deficiency at the hi^er values of the variable 
are quite possible. 

1.31 Straight infinite duct . - In order to satisfy the equation of 
motion in a bend followed by a straight duct we will match the solu- 
tion given by equation (1.26) with 


00 

i(o)t+S x+C y) 
n n’' 
e 

n=-oo 



(1.31) 


which represents an aggregate of plane waves in a strai^t tube , x 
being measured along the axis of the tube, y is perpendicular to 
it and y = 0 is at the center of the tube. The coefficients £ 

n 

and ? are to be determined, 
n 

The Helmholtz equation in Cartesian coordinates reads 





(1.32) 


Differentiating equation (1.35 and substituting into the 
Helmholtz equation we obtain 


- = 0 


£2 = k2 - 
n n 


or 
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Rewriting equation (1.31) 


<l> 


00 \ 

\ ' / +ix^ -ix i(? y+iot) 

V%=2_03 ^ 


n=-“ 

or more generally 




iwt 
e \ e 


ix 


n=-oo 


1? y -1? y 

n n 

P er + Q e 
n u 


+ e 


-ix 

' n 


p. . n 

e + e 


(1.33) 


Equation (1.33) represents waves going in the positive and the nega- 
tive directions of x and- y. Differentiating with respect to y 
and applying, at Ximits y = ±L, the boundary condition of perfectly 
rigid wall, we have 


±Z -i? R 

Pe“-Qe 
n n 


-i^ L +ic L 

Pe “-Qe ^^=0 
n n 


and similar homogeneous system for and Q^. 


The consistency condition is 

2iL -2iL „ 
e - e =0 


which gives 


sin 2c L = 0 
n 


T* w KLTT 

2 ? L = niT C = — 

n n 2L 


n = 0, 1, 2, 3, . V • 


(1^34) 
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Noting that 


k = 


u _ 2ir , 
c ~ X 


r = 21 > _2L > 2£^ , 

2L - 2L X ^ 

= 1» 2, 3 (1.35) 

X being assimed much greater than L. 

Substituting the values of P^, and ^ , from the com- 

patibility eq.uations into equation (l.33) we have 


'<j> = e 


-cot 


CO 

E 

n=0 


cos 


ijj^Cy + L) 


IX 




E' e 
n 


+ E" e 
n 


-ix 




and, writing out explicitly the solution for n=0 



If we now assume the straight tube to be infinite we must elim- 

i 

inate the exponentials with positive powers of x. With = E^ = 0 
we obtain . 


(j) = 




cos 


+ l) 



(1.36) 


/ 

which is the final form of 'the equation of motion in the straight 

J 

duct . The first term on the ri^t hand side of this equation 
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represents a plane wave. The second term is an infinite set of 
attenuated vihrations. Differentiating eg.uation (X.36) with respect 
to y will give expression for transverse vibrations. When dif- 
ferentiated, the first, permanent term yield no component. There- 
fore, only the set of hi^er modes produces attenuated transverse 
oscillations. This result discloses that the considered long waves 
in a strai^t duct must be plane to propagate. 

Based on the above then, a general incident long wave will pro- 
duce a traveling wave and attenuated transient components. This 
condition will remain until the initial wave straightens out and the 
transverse components diminish to zero. The initial wave then 
travels as a plane wave of reduced amplitude. 

1.32 9^ radians bend followed by a straight infinite duct.- To anai- 
lyze. the motion in the system bend-straight duct, it is necessary to 
determine the (in the equation of motion in the bend) and the 

(in the equation of motion in the straight tube) using relations 


1 3 d> 3ij)- 
r 36 "" 3x 


and 


It _ 

3r 3y 



(1.37) 


that is , the continuity conditions at the junction bend-straight 
tube . The velocity components , with 


L = 


Kg - 




2 


and y = r - 


2 


are : 




iiMultiplying equation (1.38) by and using the orthogonality properties established earlier 


we obtain 




(1.40) 


The integral in the denominator is given in the appendix 3. To calculate we multiply 
eq,uation (1.39) by sin[?^(r - R^)J and by orthogonality of the sine functions 



1 


(1.41) 




( 1 . 43 ) 
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al / J 

J r 


The Integral / dr is not directly obtainable. Let ns- consider 

the Bessel equation 


F" + ^F' + - ^)f = 0 


Multiplying by r and integrating (-with rF' 


R, 




» 0 , 


k' 



R^ /»R. 


rF dr = 

Ej_ R^ 


m 


dr 


/ 


An integral / rF ’ dr can be e-valuated by a general fonaiila 

m 
o 


(t)db = X r:^— 7 — ■ . - Re(p +v) > -1 
V ^ K! (v + k + 1) (y + V + 2k + 1) 


The integral 



R. 


F^ sin 
m 


r 




dr 


R, 


vill reduce , by integration by parts , to the integral 





R. 


P cos 
v_ 


m 




dr 


(1.44 ) 


R, 


which is already present in the equation (l,40)» 
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The evaluation of integral 



(1.-45) 


is relegated to appendix 4. 

The form of the solutions and of the expressions for the con- 
stants obtained so far are too complex to allow detail conclusions. 
However, equation (l.30) which defines the potential in an infinite 
bend indicates that there will be both a steady field of transmitted 

I 

velocities and a decaying field. Furthermore there will be a field 
of standing radial vibrations throughout the bend. Similarly, eq.ua- 
tion (1.36) for the straight duct shows that there will be both the 
transmitted and the decaying velocity fields in the duct and that 
the radial oscillations exist only at the inlet. You will note that 
these equations do not permit the evaluation of the amplitude and 
distribution of the oscillations. However, expanding the functions 
and integrations will simplify the algebraic expressions and result 
in numerical data. The alternating series, described by the summa- 
tion terms, converge very rapidly and will yield accurate data that 
will allow a physical interpretation of results. 
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2.0 INTEGRATED SOLUTION 

Before any numerical application of the derived equation is 
considered it is worthwhile to analyze the functions present in 
the final expression for the velocity potential 






lut 

sin(iTv ) 

yn 


C V 
m m 


sin V B + D 
m m 


cos 


V 0 

m 




;j^ (kR^) 

_ m "" m 


(kr)J^,(kR^) (1.26') 

' m m J 

The discussion will hegin hy an analysis of the Bessel functions, 
followed by an analysis of 

F (r) = J (kr)J' (kRj - J (kr)j' (kRj 

V V -V 1 -V V 1 

m m in mm 


evaluation of coefficients C , D and coefficients E for the 

m m n 

two typical duct systems to be concluded by interpretation of the 
final equations of motion. 



2.1 The Characteristic Functions F and Coefficients C 

— " , — . - y jjj — 

m 

2.10 The real root solution. - Let us consider the solution of F 

— y 

o 

in equation (l.26) pertaining to the real root v^. Using the series 

expansions for (kr) » J (kr) and for (kr) and J' (kr) as 
o o o o 

given hy equations (’l.l2,) we obtain: 


F . (r,R, ) = 


sin(TTV ) 


k(kE^)v^Tr(l - y^) ‘ 


V 

(^) * ‘*'’0 - 

V 

* (^) <‘'>^>''('’0 - '■o * 

* (^) " ‘'’P 


(1^) ° <"'*'>^('’0 + "I - 


In an direct application of this equations we would neglect 
terms containing the second powers of (kr) and (kR^^) which are small 
of the fourth order. However, neglecting terms in (kr) would re- 
strict the characteristic function to constants and terms with log- 
arithms of r. Consequently, to maintain generality we will keep 
the terms in (kr) during analysis. 
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Using power series for the exponential functions 



and neglecting terms in and terms containing products (kr)*^ 

2 

(kE^) which .are of the same order, we obtain 


o 


sin(iTV ) 
0 

4Tr(kR ) (1 - v^) 
1 o 


X 


j^j^-8 4 8v2 - in ^ I 2vl(kil^)2 in 



This general solution of the Bessel eq.uation well werifies the 
differential equation in a Tery wide range of radii ratios. For 
a = ^2^^1 1 to k the error is negligible. For a = 10 the 

error is approximately 1 percent. 

This equation, if greatly simplified hy elimination of products 
of small terms becomes 


F 

V 

o 


(r,R^) = 


sin(TTV ) 
o 

kir(kR^) 


- 8 + 8v^ - 2(kR, 
o 1 




+ 2(kr)^ - k(kR^)^ 


(2.3) 
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Equation (2.3) » when applied to the differential equation, still 
yields a very satisfactory expression for at r = Eg namely 


v2 


- 1 


In a 




-2 


= (kE^) 


2 - 1 
2 In a 


Eor numerical calculations F is, -with hi^ accuracy, equal to 


2 sin(nv ) 

ir(fcE ) 
o 1' 

To obtain C^, we have, by equations (1.24) and (1.25), 
the expanded expression for F 

o 

V 

C 2ii(kHj(l - v^) X 

O V 1 0 

o 


(2.4) 
and using 


{[‘' - '’“V 


) + 

O 


(a - 1) 




2(kE^)' 


4v2 

o 


a In a + 2 v 2 a(ln a)^ 
o 


^6 + 8(kR^)^ In a + 4(kR^)^ (a^ - l) + 8(kB^)^ (in a.f + f(kR^)^ • 


■which, if the small of the second and of the hi^er orders and neg- 
lected simplifies to 


w{kR^) 




a " 1 
In a 


(2.5) 


o 

2.11 The imaginary roots solutions . - Let us consider the Bessel 

functions of pure imaginary order in the equation (1.26). Use will 

be made of expressions for the Bessel functions with real argument 

and of imaginary order and tables of functions which have been pub- 

( 20 ) 


lished by Buckens 
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The very basic relations are 

F^(kr) 4- iG^(kr) = 2^^r(l + iy)J^^(kr) (2.6) 

where iy = v and 

G^(kr) = A^(kr) sin [y In (kr)] + B^(kr) cos [y In (kr)] 

F^(kr) = A^(kr) cos [y In (kr)] - B^(kr) sin [y In (kr)] 

The functions A (kr) and B (kr) are A (kr) = 
y y y 

OO ^ 

®y ~ g = 0, 1, 2, 3. . . uid 

(?)^ = 1» (s)^ = 0 and 

r - ^j2g-2 ~ ^^2g-2 

^2g k5(y^ + g^) ^2g ''4g(y^ + g-^) 


Neglecting the fourth and higher powers of (kr) we have 


'^(kr) = eostu In(kr)] - 4 ' ( 2 ^V ' ^ ' ^Y y In(kr)] 


G (kr) 

y 


_ (kr)^ 

~ ^ - 4(1 + y^)_ 


sin[y In(kr)] + 4(2.^+"y^ ' ) ' In(kr)] 

(2.7) 

Substituting into equation (2.6) we obtain expressions for 

J, (kr) and J . (kr) for (kr) << 1. In an abbreviated notation 
ly -xy 

they are 

F^(kr) + iG^(kr) = 2^^r(l + iy)J^^(kr) 


F (kr) - iG.(kr) = 2“^*^r(l - iy)j . (kr) 
y y -xp 
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The left hand sides of the above equations are complex conjugates . 
Since 

r(l + iy) = r(l + iy) 

"we conclude that ^^(kr) is complex conjugate of 2 ^^J^^(kr) . 

2 

In equations (2.T), for (kr) « 1, terms in (kr) may be .neglected. 
Expressions for (kr) and J ^ (kr) simplify as follows 

cos [p In (kr)] + i sin [y In (kr)] = 2^^r(l + iy)J^^(kr) 

(2.8a) 

cos [y In (kr)] - i sin [y In (kr)] = 2^^r(l - iy)J_^^(kr) 

(2.8b) 

a 

Let us substitute into equation (i.8a) the value y = m ^ 

and consider the term sin [y In (kr) ] for kr — kE^^ and kr = kE 2 

We obtain 



The second expression is dephased by mir with respect to the first 
one. We conclude that for m = 1, [y = Tr/(ln a)], J^^(kE 2 ) is the 
negative of because it is reflection of 

origin. For m = 2[y = 2ir/(ln a)], J^^(kR 2 ) will be niimerically 
equal to J^^(kR^) but will be rotated hy 2 tt. The phase angle of 
the solutions for (kR^-) and (kRg) hy even and odd numbers of ir im- 
plies that in the range of radii from R^ to Rg there will be a 
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number, increasing with m. of pure imaginary or simply real 


J. (kr). 
ly 

Taking the deriwative of the expression (2.6) with respect to 
the argument (kr) of the Bessel function we obtain 


xy 


2^^r(l + iy)jj„(kr) = sin[(y In (kr ) ] |2?^(kr ) + y|^(kr) 


TkFTj 


+ .cos [y In (kr)] [y?g(kr) - 252(kr)] 


■(- 


'[■ 


+ i ( - cos [y In (kr) ] j 2C2(kr) + yCg^^^) “ 




+ sin [y In (kr) ] |y?2(kr) - 2Cg(kr) 


} 


All terms containing ?2(kr) and g^iVir) are small when com- 
pared with y/(kr) and in the first approximation may be neglected- 
The res\ilting eq,uations are 


2^^r(l + iy)J^^(kr) = - ^ (hr)] + i cos [y In (kr)] 

(2 .9a) 

2~^^r(l - iy)JV^(kr) = - sin [y In (kr)] - i cos [y ln(krl] 

(2.9t) 

By considerations identical to considerations applied to expressions 

(2.8) we conclude that the derivatives' with respect to argument of 

the Bessel factions of pure imaginary order J^^^(kr) and J|^(kr), 

for (kr) « 1, are complex conjugates. 

Using equations (2.8) and (2.9) we form for m 5^ 0 

m 
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F 

V 

m 


y sin[y In(kr) ]sin[y In(kE^)] + cos[y In(kr)] cos[y In(kK)] 
2 ^ m m m X la ml 

r(i + ±)ij r(i - xy„) 


m 


m 


( 2 . 10 ) 


Since 


r(l + iy) r(l - iy) = ^ 


TT]i 


sinh (iry) 

and using trigonometric transformation we have 


2 sinh(Try ) 

-r, _ . m 

% “ Tr(kR, ) 

m 1' 


y In 
m 




( 2 . 11 ) 


Equation (2oll) shows that F is a pure imaginary number and is 


m 


a function of radius r. To check the degree of approximation used 

* 2 

when neglecting -terms in (kr) we retrace our calculation and estab- 
lish F (r,R^) with better accuracy. It is 
m 


sinh(iTy ) 


F 


= 1 


m 


m 


4(1 + y^) Try (kR) 


X 


m 


m 


{■ 


2y 


m 


4(1 + y2) + (kE^) + (kr) 


4(kE^ )^ + 2y^(kR)^ - 2y2(kr) 
1 m m 


'] ““(-m I t) 

'] i^)] 


Except for the first term in the expression of the cosine of the 
logarithm, all other terms are very small and are rapidly decreasing 
with increasing m. The simplified expression obtained by elimina-' 
tion of terms (kr) , which are small of the fourth order, retains 
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the character of the exact expression and may he considered as a 
very good approximation of the exact equation. 

To calculate the constant pertaining to the imaginary 


roots we evaluate the two integrals, | F dr and / — 

m J 


F^ dr, 
V 

m 


using equation (2.10) for F . This expression contains the vari- 

m 

able r in functions sin[p In (kr)] and cos[v In (kr)]. Before 
proceeding, we establish the necessary integrals starting with the 
well known integral 


A 


X X 

sin(ln r)dr = ^ r sin(ln r) - -g r cos(ln r) 


They are 


/ r^. sin()i In r)dr = =r (p-+ sin(v In r) 

J [(p + 1)2 + y2j L 


(P+I) / T \ 

y r cos(y In r) 


/ 


r^ cos(y In r)dr = 


|jp + 


1 Tr 

1)2 + y2j L 


sin(y in r) 


+ (p + l)r^®"^^^ cos(y In r) 
where p is unrestricted. 


sin^(y In r) , 1 , 1 . fn \ 

L dr = — In r - t — sin(2y In r) 

r 2 4y 


f- 

cos2(y^ln r) = |. jn r + sin(2y In r) 

J sin(y In r) cos(y In r) eir = - cos(2y In r) 


( 2 . 12 ) 


(2.13) 

(2.14) 

(2.15) 



[p sin(2p In r) + cos(2p In r)] 


(2.16) 


Up[l + p^] 


j r co=(, 1» r) sln(l. In r)dn = - 
y^rW(n In r)dr - i f °id(2n In r) ^ ^rj2 _ 


(2.17) 


The maxn lines of calculation of are quite sample but lengthy. Omitting the algebraic intermediate 

steps, the result is 


v^Ti R^(kR^) 

^ )Tii sinh(irp f ^ 
m m 


^ 2(kRi)"^_ 2(kR ) 

v“6h (a cos mir - l) (2 + p + p^~)[a cos mir - l] + |p3 _ 3 y II cos mn - ll 

1 + p2 ^ ' 32 + yzLm mJL J 


m 


m 


J 


1 

1 1 

^1 ^ 

V 

1 + 

' 2 + — y 
p‘‘ 

m m 

1 + 



Jj 


and retaining the large terms only 


C V 
m m 


= i 


\ (a cos(mTr) - l) 

sinh(p ¥) In a(l + p^) 
m 


( 2 . 18 ) 


and 

Q _ ^o^^l^^l (a cos (mir) - l) 

m m sinh(p ir) (1 + p*^) 
m m 


(2.19) 



2 . 2 The, I n fin ite . Bend. 


2 , 20 Equations ♦ - The integrated equation for the velocity potential <j>(r,0,t) fo'r an infinite 
bend may now be written directly by substitution of the derived expressions for C^, and 
into equation (1.30). 


o m 

The velocity potential is 


(|) = - i 


V Rt> 
o 1 

2v 


l(,(0t-V 0) 
e 


a - 1 
In a 


-2 + (kr)^/2 + 2v^ - 2(kE^)^ 


- ^ 1“ - I) 

J. « 


+ 


V R, 

o 1 

4tt 



a cos mir,- ,. ! /g^ + ^ 2 ) + (^r )2 

my (1 + y2)? ] “L “ ^ 

m m V. 


+ (kr) 



In 



+ 


4(kE^)^ + 2y^(kR^)^ 


2y^(kr) 



and without terms small of the second order 


(2.20) 


((i(r,0,t) 


i 


V . 

“R 

V 1 In a 
o 


i(o)t-y e) 

O 

e 


0^1 icot \ 

— ^ z_ 


a coB(mTr) - 1 
m(l + y^) 




( 2 . 21 );' 


<-n 

hO 



The tangential vibrational velocity is 


, ^ , iwt 1 a - 

v(r,e,t) = -e 


iv er 

e “ [- 


2 

2 + (kr)^/2 + 2v| - - v|^ln 


l2\ i«tV 

' “ 4iTr ^ / 


m=l 


m a cos miT - 1 

m(l + v^f 




[4(kE^)^ * 2n2(kr)^jsik^p_^ In |-^\ 


(2.22) 


The radial vibrational velocity is 


V 3.(a)t-v 0) 

/ ^ \ . o 1 o a 

u(r,e,t) = -1 -^r- — e 


2v 


V ~ - l k^r - (kE, )Vr - 2v^ In |- / r 
lna[_ 1/ o \/„ 


+ e 


ia»t 


00 

V R, 
_o__l ) 

Uv r / 


( a cos(rav) - l) 


m=l 


m.(l + p^)‘ 




U (1 + y2) + (kR )2 + (3j3,)2 

m. X 


- ity^(kr)^sin^y^ In + |^(kR^ 


+ 2 y 2 (kR _ 2y2(kr)^ - 4(kr)^ 
ml m 


cosly^ In - 


Ln 

LO 


(2.23) 
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The real parts of these equations fully describe the propagation 
of acoustical waves in an infinite bend. To better visualize the 
dependence on radius of the final equations of motion, figure 2 
shows cos ln(r)/(Ej^) which appears in equations (^.22) and 
(2.23). This function is not very common, ¥e note that the traced 
curves • are shifted off-center due to the influence of the log- 
arithm in the argument of the cosinus. 

2.21 Discussion . - In the expression for the velocity potential 
the term which rapidly decreases, with increasing 0, is ir/2 out 
of phase with respect to the other term' which is a true wave of 
amplitude 


V _ 

o a - 1 

2v 1 In a 
o 


:-2 + (kr)^/2 + 2v^ - 2(kR, f' 
\ o 1 




In 


r 




Since kr << 1, this amplitude is, in a first approximation, in- 
dependent of radius and the progressive wave is nearly plane. The 
tangential vibrational velocity has both the permanent wave term and 
the attenuated vibration terms in phase. The permanent wave term 
is basically inversely proportional to the radius while the atten- 
uated, disappearing term has a more involved radius dependence. 

The radial vibrational velocities also have a permanent and 
an attenuated term. The permanent stationary vibrations are ir/2 
out of phase with the tangential components. The attenuated oscill- 
ations are in phase with the tangential velocities. 




Figure 2. -Cos[ti,^ ln(r/R^)] for m = 1, 2, 3, and 4, 
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A vectorial addition of the radial and tangential velocities 
yields motion's vihrational velocity and the direction of the wave 
front. 


The propagation of waves in a curved duct is thus profoundly 
influenced by the curvature of the conduit. The amplitude of 
vibrational velocities is a function of both inner radius and 

of the radii ratio a = The tangential vibrational velo- 

cities are almost exactly inversely proportional to radius, that is, 
velocities follow the distribution of a potential vortex. The 
radial vibrational velocities are approximately two orders of magni- 
tude smaller than the tangential velocities . Their radial distribu- 
tion is governed by the three terms in the square bracket of the 
equation (2.23). The expressions derived for vibrational velocities 
in a, bend permit a thorough discussion of the motion. The results 
obtained will be discussed with the aid of figures 3 to 8. 

The tangential vibrational velocities of the propagating wave, 
as described by equation ^.22), very closely approximate a hyperbolic 


distribution across the channel's width. Figure 3 discloses that the 
vibrational velocities for a = E^R^ up to 2.5 range from values 
close to 1.5v^to approximately 0.5"v^* It is noted that value of 
Y = always falls on radii smaller than ducts centerline radius 
R^ = (R^ + 'B.^)l2. For example, for a = 2, v = v^ is at 
r^^ = (a - l)/ln a = l,kk2'J while Ern^R^ = 1.5- Performing 

the calcvilation for different values of the a, we obtain the re- 


lationship r /Rm as a function of the parameter a 

Y=V^/ 






Figure 5 - Attenuated tangential vibrational velocities for three angular positions in a bend 


Phase velocity ratio, |Sj/x 
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This is shown in figure 4. The radius on which v = decreases 
proportionately to the increase in the sharpness of the bend. 

The phase velocity of waves is also affected hy bending of 
ducts. For a bend, the phase velocity is S = obtained 

from (uit - while for a straight duct it is x = to/k, as given 

by {fjut - kxj . To compare the two velocities we average the tan- 
gential phase velocity, 0r = s over the duct width and obtain 

M kR 
= — — R = — ^ 

' ’ V m • e V m V 

o X 0 o 

This ratio is traced on figure 4. Clearly, the phase velocity in 
bends is always higher than in straight ducts . Furthermore it in- 
creases with increasing duct sharpness as measured by a = Rg/R^. 
Therefore, if equal passage time in a straight and in a bend duct 
or phase equality are desired, the bend pipe must be made longer 
than the straight pipe. 

The attenuated tangential vibrations which characterize change 
from motion of plane waves in straight duct to motion in a curved 
duct are examined in the next two figures . Figure 5 gives results 
of a sample calculation illustrating the behavior of those oscilla- 
tions for a duct of radii ratio a = 2. The vibrations are basic- 
ally of low amplitude. Even close to the piston, at 6 = they 

are one order of magnitude smaller than vibrations of the non-damped, 
propagating wave. The radial distribution of these oscillations 
changes significantly with wave angular position in the duct. At 
e = Tr/4 these oscillations are reduced to a very low level and are 
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nearly laniform across the width of the duct. Figure 6 shows the 
same oscillations calculated for three different duct widths hut 
with a single angular position of 0 = tt/ 4. The curve for a = 2 
was taken directly frim figure 5 for comparison purposes . The two 
other curves indicate that the decaying oscillations are much more 
pronounced and extend farther when induced in wider ducts. 

The radial .vibrational velocities , which characterize motion 
of waves in curved ducts, for long waves, are of low amplitude. The 
permanent, standing oscillations are shown in figure T. They are 
calculated for duct radii ratios of 2, 3 and k. Generally the amp- 
litude of these oscillations is low, approximately two orders of 
magnitude smaller than the tangential velocities. The radial dis- 
tribution is characterized by the lack of symmetry. The maxima of 
curves are shifted toward bend's inner wall. 

This phenomenon is even more pronounced in the case of the non- 
propagating, damped, radial oscillations at the curved duct inlet. 

The amplitude and radial distribution of these oscillations is shown 
on figure 8 for three values of parameter a and at 6 = ir/k. For 
a = -2 the permanent and the vanishing oscillations are of approxi- 
mately the same amplitude. For a = 3 and U the decaying oscilla- 
tions are about twice as large as the radial oscillation of the prop- 
agating wave. Therefore, the process of decay is much slower when 
ducts are wide. 

The initial conditions , that is , the distribution of velocities 
on the vibrating piston, determined the type of the obtained solution. 




Figure 6. - Attenuated tangential vibrational velocities for three bends at 6= tt/4 



Figure 7. - Standing radial vibrational velocities Figure 8. - Attenuated radial vibra 

in bends for three Rj and three values of a tional velocities at bend's inlet 


6i 


To verify the influence of the initial conditions let us compare the 
simple initial distributions in which the vibrational velocities at 


the piston vill be independent, proportional and inversely propor- 
tional to the radiiis. Calculating the integral | v F (r)dr of the 

J m 

nominator of the constant Cm we obtain the following results 
Initial Amplitude of the Amplitude of the "m"^^ Decaying Wave 


Velocity Progressing Wave 
Distri- is Proportional 
but ion to 

V f + 0(kr)^ 

o In a ' 


is Proportional to 


(a cos mir - l) + 0(kr)‘ 


— + 0(kr)^ 2 — t . (a^,.cos mtr - l) +/0(kr)^ 

^ m^ir^ + 4 (in a)^ 


zero + 0(kr)^ 


The first possibility, that is the particle vibrational velocity at 
the piston are independent of the radius , has been already discussed. 
The second possibility resTilts in the amplitude of the wave and of 
the decaying vibrations strongly increasing with channel radii ratio 
a. The third possibility (potential vortex at the inlet) results in 
a wave of almost constant amplitude , independent of the parameter a 
and with vanishingly small damped oscillations at the inlet. The 
assumption that the tangential velocity at 6=0 be inversely pro- 
portional to the radius results in a remarkable simplification of the 
equation of motion because virtually only the imdamped simple pro- 
gressing and weak radial waves are present. Furthermore, the pro- 
gressing waves are essentially independent of bend's radii ratio. 
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The physical meaning of the derived equations may now he dis- 
cussed. We will make a direct parallel between solution of the 
Laplace equation = 0 for inviscid, incompressible fluid and 

the Helmholz equation V^<() + k^<j) = 0 for inviscid , compressible 
fluid. One of the solutions of the Laplace equation, written in 
cylindrical coordinates, gives the potential vortex. It is charac- 
terized by absence of radial velocities and, at every 0, the tan- 
gential velocity distribution is the same. 

In the case of the Helmholz equation with v^ = const, or 
proportional to the radius, equations indicate that long waves only 
gradually adjust themselves to the curvilinear boundaries vintil a 
permanent pattern of a free vortex type (in tangential vibrational 
velocities) is obtained on which vanishingly small components of 
linear and logarithmic nature are superimposed, and the previously 
discussed set of small radial oscillations is formed. 

We conclude by saying that solution of the Helmholz equation 
for a curved, two dimensional duct results almost exactly in a 
potential vortex in- the tangential velocities with superimposed 
pattern of .small standing radial oscillations. The presence of 
these radial vibrations distinguishes the Helmholz equation from the 
Laplace equation type of motion. 

2„3 9^ Radians Bend Followed by a Straight Infinite Duct 
2.30 Motion in the bend. - To calculate the coefficients D which 

ja 

determine amplitude of vibrations near the outlet of a bend,, we will 
proceed with integration of expressions established in the general 
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solution. Equations (l.40), (l.4l) and (lAs) indicate that motion 
in a hend depends on the parameters of the strai^t duct that follows 
the bend. Obviously, motion in each segment of the ducting depends 
on motion in the entire system. The derived relationships are quite 
involved and simultaneous solution of these equations results in an 
infinite matrix in m and n. The coefficients depend on co- 

efficients. C , E and E . However, by substitutions, it will be 
m’ o n ’ 

possible to express D in terms of known C only. 

m m 

To simplify the writing the following symbols will be .used. 

For the integral of equation (1.45) 



?„(r - B^) 


dr = m ^ 0 

Emn 


(2.24) 


The coefficient E , as given by equation (l.43) is 




m2 iC 


dr 


The first term, which will be designated by S, consists of two 
parts. One large and independent of the parameters of the bend and 
a series of very small terms which are functions of and a. 

The second term consists of a sum of very small teims only because 
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the main term of the integral 



(r)/r dr 


is 



In an expanded form and with the help of simplifying symbols , equa- 


tion (2.25) reads; 



Icoth^y^e^] 

t - tanh(iJ^e2) 

/ V 

+ X 2 

coth^ygeg^ 

) - tanh(y2®o) 


+ + XgPg * . . . 


where 3 = 

m 


D 

m 


i C 


+ tanh 


m 


(V2) 


(2.26) 

(2.27) 


We note- that for than apprximately 4, coth 4 ~ 

tanh 4 ~ 1. Typically for 0^ = ir/S and a < 5 this approxima- 
tion is V63ry satisfactory. In alL practical applications of long 
wave motion in a 90° bend a simplified equation 
V -iv 0 

E =-i ^ e ° (2.28) 

O 


may be used with hi^ accuracy. Equation (2.28) results from equation 
(2.26) when small terms of the second and higher orders are neglected. 
By equations (l.40) and (l.4l) a single expression for D 

I 

may be obtained. It will be convenient to establish the ratio 




(2.29)' 

/ 

where xs given by eq.uation (2.26) and where 'D is the given particxilar D taken into con- 
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sideration. Equation (2.29) cannot be simplified any further. It contains an infinite matrix in 
m and n because each depends on all and on the parameters of the strai^t duct. Ex- 

panding the matrix for m = 1, 2, . . . and for n = 1, 2, . . . and introducing the abbrevi- 
ating symbols given by equations (2.26) and ('2.27) we obtain 







-1 


sinhl u_ 6 _ )v , 

Cv dr 

\*^1 2 } 0 J 



_ k^'" A . - ~ 
I 1 V 1 Fli-TT 


iC 


sia(.vj ' (>' 1 ® 2 ) '^: 


’Fll 


iC, 


+ TT 


2 sin(iTV^) '^°®^l''2“2/‘F21 


^(^^2^2) ^2 


+ • . 


+ ?2 V *^2 ■ h 


iC^ 

^1 sin(Trv^) ^F12 



- 00 th 



« • 


^*^2 / \ 
- ^2 sin ' (^vj 


F22 



( 2 . 30 ) 


Grouping, rearranging and adopting additional symbols to designate the constants of- the" 


matrix, we have 



-i- coth ^2.^2 


00 

'L 


Q , 


^ " ^1 Z-^ ■'In 

n=l 


- X 


■11 






B, + coth ( 

li e ) 

- tanh I 


1 \ 

ri 2/ 

ri 2/ 

k “*** ( 

.^2^2) 

- tanh 1 

h^2[ 


+ ,+ 123^62 ■" " 31^3 * 


(2.31) 


Bg + coth (^2^2) ~ (^2^2] 


Ag - Zg 


CO 

-z 


Q , 


'2n 


n=l 


CT\ 

03 


~ Xn 


- Xr 


B^ + coth ^U2®2’) " tanh^U^Bg^ 
Bj.H- 00th (V282) - t“h(B2®2) 


+ nj^2*l ^'"22^2 "32®3 * ' ' • 


Where A^ designates the second term in each row of equation (2.30)j Z^ are used for products 


ik 


Z ~ — "q — and Qi stand for terms in the third line of equations (2.30) . The x. 

0 ' ffl 2 


mn 



are mTJlt~ 


are multiplied by the factor in m in front of each square bracket and the 

ipliers of all B in equation (2,30). Finally, by rearrangement, the matrix receives a symmetrical 
m 

form. 



and similar expressions for the other rows of the matrix. The derived solution is very general. 
It applies to radii ratios of the bend in the range from 1 to 10; it has no limitations as to 
the length of the bend which may range from a small fraction of tt to an infinite coil. This 
wide range of applicability was obtained by retaining in the solution, terms which are two 
orders of magnitude smaller than the principal terms . Now, reasonable limitations to this range 
of bend parameters will greatly simplify the final matrix. For bends with Sg greater than tt /4 
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or tt/2 and with a <2.2 or a<5 respectively, and neglecting 
small of the second order, the expressions in the square brackets 
drop ’OUt, Furthfermd're, the 'Z become - Z and A , 0 and x 
tend to.be negligibly small. ^ .The left aide of -each row of the 
matri'ix' b'e’come -a con’s fant -throughout ' the matrix and the-‘right side 
further '-simplxfies the matrix, equation '(2. 32a)V^reduce^ to 


"2 “ K ^1^1 ^^21^2 “ ’^31^3 

-2 == 4- (n*2 4- i)3' - ' 


(2.32b) 


where "prime" on all variables indicate that they have been simpli- 
fied by elimination of small of the second and higher orders. 3' 

m 

are now defined by 



= - 1 - 3 ’ 
m 


-iv 05 

O d. 


This simplified form of the matrix which determines constants of 
motion in a bend is quite satisfactory for a wide range of applica- 


tions . 


D 

With the expression for ^ ^ " already established and using 

m 

the adopted abbreviating symbols, equation (l.4l) for is 


simply 


E 

P 


2 "o ? 

p-ir Ysin( Trv ) 




00 


cosh 


% 

(l'm®2) 


m=l 


m sintrv ) 

m 


F’ (r)sin 
'^o 





(2,33) 
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where is the given particular talcen into consideration. 

Kow equations for the velocity potential and the vibrational 
velocities are fully determined. However, expressions for <j>, v 
and u respectively, cannot be written in a simple form because of 
use of series expansions of the Bessel functions and the involved 
expressions obtained for constants . 

However, should we consider ‘large terms only [neglecting small 
of the order of (kR^) ] and bends larger than it/4 so that 


1 ^^2 

sinh wSg “ cosh “ — e 

we may form simple equations which with satisfactory accuracy des- 
cribe the motion in curved ducts . In fact , these simplified equa- 
tions make it possible to discuss several features of these rela- 
tionships which would be otherwise obscured by large number of 
small terms which add to the exactness of the solution but not to 
the determination of the essential features of propagation except 
for the presence of the permanent radial standing vibrations. These 
vibrations, which characterize the motion of waves in bends, are 
described by terms small of the second order. Since these radial 
vibrations have been adequately discussed in connection with equa- 
tion (2,23), there is no need to dwell on this subjject again to any 
length. 

The final approximate expressions for propagation of long waves 
in curved ducts will be based on equation (l,27) which may be written 



as follows 


iC F 
m V 

1 

2 sin(7rv ) 
m 

Introducing tlie expressions for 5 F^ and for [ j gi'ven respectively "by equations 

^o m V m / 

( 2 . 11 ) and (2.32a), we obtaxn 


iC F 

o V 


i(wt“V 0) 


~ sin(Trv J 


..-I 


m=l 



^(r,9,t) = i 


. a - a , l«t 


V In a 
o 


+ 6 


0^1 a C' 

—L ^ 


m=l 


cos(nnr) - 1 
m. 




cos ly„ In — I X 


-y 6 

-2e ™ + B e 

m 


-V r 

o 2 \| 


ye -y 9 ^ 

e , + e ; 

y 0/a 
2 


(2.34) 


-iv 0 p 

where the term e ° results from evaluation of E at x — 0 which is at 0 — 6p and 

n ^ 

not at 0 = 0 and could he avoided should the coordinate 0 of the bend he counted clockwise from 
the junction hend-straight duct. 

Differentiating with respect to 0 we obtain the tangential vibrational velocities 
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(2.35) 


The first two terms of equation (2.35) are those of equation (2.22) of the wave motion in an infinite 

hend. The first term is a simple wave, the second the attenuated vibrations near the inlet section 

of the bend. The third term represents attenuated vibrations in the negative direction of 0, an 

adjustment of wave motion to discontinuity at the junction bend-straight duct. This term is identi- 

-iVo02 

cally zero at 0=0 and increases to a maximim’' at 0 = 02» Because of the presence of e 
these attenuated vibrations are in phase at 0^ but out of phase of motion at 0 < 0g. Actually, 
this phase difference represents just time shift. We have ‘^qQ 2 ~ "*^2 ■tiine required 

for the wave to reach section Equation of motion in the straight duct have time counted from 

t = 0 at X = 0 or 0 = ©g while actually the time count should have been started at a later 

2’ 


time t = t, 
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The radial vibration velocities are obtained -by. differentiation 
with respect to radius 


u(r,0,t) 


i(ot 
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(2.36) 


This expression for radial vibrational velocities do not have any 
permanent wave term which was given by equation (2.23) when small of 
the second order were considered. The damped vibrations described 
by equation (2,36) exist near the inlet and near the outlet sections 
of the bend only although, theoretically, they extend over. the whole 
propagation space. The amplitude of the inlet induced velocities 
is different than an^jlitude of exit Induced velocities (factor 1 vs. 
6/2). Also the velocities near the inlet section are of opposite 
sign than velocities near the exit. 

'The motion in the straight duct may now be analysed. The co- 
efficients 3 are given by equation (2.33>. Substituting the 
^ D 

established approximate values for and and simplifying 

m 


we obtain 


E 


R, 


o irL 



m=l 


a cos(mTr) - 1 . „ 

m(l + v^) mp'^m 


(2.37) 
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where are given hy series and matrix in appendix 4, and 

are numerical coefficients. Substituting equations (2.28) and (2.37) 
into equation (1.36), the velocity potential in the straight, infinite 
duet is approximately 

i(ut-v 0.) V 
j,\ o 2 . o -ikx 

4>(y,x,t) = e ^ ® 


2v E. 


■ E, C 
o 1 \ n 

ir^ Z_ n 


cos 
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?^(y + Ii)|e 
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Z A a cos(m'iT) - 1 o 
mn ~(1 + y^) 
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(2.38) 


Differentiating with respect to x and simplifying, with 

+ k^ - ? , the axial vibrational velocities are 
V n n 
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The normal vibrational velocities are 
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-e 


09 



n=l 


n sin 


+ l) 



X 



A 


mn 
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(2.40) 




Hie double series in equations (2.59) and (2.-40) are rapidly con- 
verging as long as x is not zero. The time shift coefficient 
-iv 

e is synchronizing the propagation of vaves in the straight 

duct with wave motion in the bend. 
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3.0 SKEEGY.FLOW 


Currently, no formulation of energy density or energy flow has 
"been developed for long wave motion in "bends* Karpman' attempted 
to develop expressions for the kinetic and. the potential energies of 
standing waves in a torus. However, because of mathematical diffi- 
culties he limited himself to partial results in terns of acoustic 
energy density which is sound energy per unit volume. 

In systems with progressing waves it will he appropriate to e- 
valuate energy flow or power per unit area, called sound intensity, 
which is given , in most general terms , by 



The V designates the vector of the vibrational velocity normal to 
the wave front. The first integration, over the period of oscilla- 
tion, is designed to replace instantaneous values with time averages. 
The second integration is necessary in the case of wave parameters 
dependent on radial coordinate to obtain an average over a cross 
section. The integration of the averaged product pv over a cross 
section is not necessary for plane waves in straight ducts hut is 
necessary for plane waves in curved conduits. 

Comparing sound intensity in a bend with sound intensity in a 
straight duct, will verify the derived equations of motion. To he 
successful, this comparison must involve the propagating waves only. 
To exclude the attenuated waves it will be assumed that in equations 
(l.30) and (1.36) the series with terms corresponding to m = 1, 2, 3, 
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• • • and to n = 1, 2, 3 * * • are exactly zero. This will be 
very nearly true for sections of bends at angles close to Tr/2 or 
larger and for sections in straight ducts at length/width ratios of 
2 or larger. 

For the bend, the tangential vibrational velocity, as given by 
eq.uation (2.22) , for large 6 and when time averaged is 



R r* 

= - V — f - r— - 2 + (kr)^/2 + 2v^ - 2(kR, )2 

o r 4 In a o 1 


- In (In 

1 R^ 


This velocity, when averaged over a cross section by integration is 


V n = V 

' 'bend o 


5 (R_ 

o 1 


(3.1) 


where 6^ is a small quantity-of the second order. Equation (3.1) 
indicates that the kinetic energy of the tangential vibrations in 
the bend is less than the kinetic energy of 'the wave- generated at the 
face of the piston. Some kinetic energy was transferred from the 
progressing wave into the standing radial oscillations. 

Using the equations on page 37/' and the relation (2.25) we obtain 

n'lr 

the tangential vibrational velocity in the straight duct, at a sec- 


tion of the duct remote from the junction with the bend. 
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In order to obtain an average value across a cross section, it is 
necessary to replace by kr and integrate the above equation. 

The equality of and kr was explained in detail in section 1.3.- 

As you recall it simply indicates that motion in a bend consists of 
an infinite set of waves of different wavelength but of a single fre- 
quency. Therefore making the substitution and integrating results in 
the following expression. 


IpI =-£^v 

I bend k o 


Ij- - 6 (a. 



In a straight duct, time averaged pressure is described by the ex- 
pression 


IpI = i ^ E 

'^'straight 2 o 

It is noted that this expression parallels the expression for 
Kl gtraight * taking the proper substitution and performing the in- 


tegration, results in the following expression, 

in PJi 
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Examining this expression we note that the potential energy of the 
progressing wave in the hend, and downstream from, it, is also sub- 
ject to a small decrease. 

Finally, using the derived expressions, we can describe the 
flow of energy in the bend by the expression 


I 


“bend 


ptD 

2k 




2 


and in the straight duct by the expression 


_ 2 
■^straight 2k o 


[l - «^(a. E^) - E^. e^) 


The long waves under consideration are subject to a loss of en- 
ergy when going through a bend because some energy is -retained in 
radial vibrations. However, at low frequencies, waves pass through 
bends very easily. This is because the reduction in the rate of 
flow of energy is very small. 

A verification using the energy method of the approximate ex- 

i 

pressions, equations (2,21) and (2.38)', may be made directly. We ob- 
serve that the approximate expressions pertain to plane waves which 
are defined by 

p = (wt - n ■ r) 

where n is a unit vector perpendicular to the wave front. For a 
plane wave the sound intensity vector i is 


2 

— — ?— p — 

X = pv = pcv^n = n 
pc 


where the quantities p and v are already averaged over a period. 
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For the three expressions of the energy flow we have 
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which, as before, only requires substitution k = v/r to verify the 
identity. 

Finally, comparing the approximate expression for the straight 
duct and for the bend we have : 


for the straight duct 


i 


pc 



2 


where jv^| is averaged over a period. 

For the bend, averaging v. over cross section yields 

o 



V 

o 


which verifies the energy flow equations. This verification is 
based on approximate expressions in which small terns of the second 
and higher orders were neglected. However, the small terms describe 
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the radial oscillations that characterize the motion of waves in 
c\irved ducts. Since no consideration was given to energy which was 
transferred from the progressing wave to the standing radial oscil- 
lations j the energy flow in the strai^t duct equals the energy flow 
in the bend. This also equals the energy flow generated by the pis- 
ton. There are no experimental data available to quantitatively 
verify these results. The reason is the physical difficulty of try- 
ing to examine the (very low) range of frequencies which is the sub- 
ject of the present analysis. Usually the acoustics experimentation 
is carried out in the audible frequency range of 200 to 3000 Hz. 
Consequently only a qualitative verification is possible. This ver- 
ification must be based on the extrapolation of the available exper- 
imental information to very low frequency range. 

(25) 

In 1955 Lippert published data on characteristic transmission 
factor T = B exp ig of square and rounded 90° elbows. The ampli- 
tude B of the factor T is defined as ratio of sound pressure am- 
plitudes of the transmitted wave at the output plane to the sound 
pressure amplitudes of the incident wave at the input plane, g is 
the phase angle of the' transmitted wave with respect to the phase of 
the incident wave. As shown on figure 9 , transmissivity of a 90° 
mitered bend is high at Icn-r frequencies (_< 600 Hz) but close to 2200 
Hz approaches zero. An experiment with a rounded ‘90° elbow indi- 
cates that rounding improved the transmissivity over the entire range 
of frequencies and eliminated the cut-off condition. G?he dotted line 
indicates extrapolation of Lippert 's data. We conclude that at very 




Figure 9. - Transmission of acoustic waves through 
90° mitered and rounded elbows. Experimental data 
by Lippert'"\ 
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low frequencies (l to 10 Hz) the transmissivity of the two elhows is 
very hi^ hut not one-hnndred percent. There will he some reflec- 
tions and some transfer of energy from the progressing wave to the 
standing radial oscillations. At the limit, zero frequency, there is 
no wave motion and the potential flow of matter shoiild suffer no 
losses of energy. That in effect is the case of a potential vortex. 
The lower curves on figure 9 show the phase change of the transmit- 
ted wave. The square mitered elbow and the rounded elhow produce 
similar phase changes , with virtually merging slopes at lower fre- 
quencies. Again, using extrapolation, we conclude that at low fre- 
quencies there will he some phase shift. 

On the basis of these dat^ it may he concluded that results of 
this analytical study are qualitatively verified. The long waves 
penetrate and pass through hends easily. However there is some en- 
ergy being stored in the hends which comes from the- main wave. This 
stored energy is a very small, hut not negligible. 
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4.0 NUMERICAL APPLICATION 

The purpose of the numerical application is to obtain a solution 
for a conventional industrial setup which will give typical wave 
propagation pattern. At the same time it is intended to analyze the 
Bessel Pimctions of very small order and Bessel functions of pure 
imaginary order and real argument. 

4.1 Selection of Parameters 

In order to simplify the numerical calculations and to conform 
to the req.uirment that parameters must approximate a typical in- 
dustrial installation we select the logarithm of the radii ratio. 

In a, such that defined by, 

V = iw = im • ^ m =5 1, 2, 3 be an integer 

m m In a 

tt 

Selecting In a = O.6283I we obtain = 5> v = im5 and 

Xn Qd m 

a = ~ 1.8T^U. This radii ratio is reasonable and approximates 

an industrial ducting. 

To obtain the argument of the Bessel functions which will 

be a small, finite number, and a wave constant k and radius 
typical for industrial application we select 
k - 0.1 m E^ = 0.2 m 

It follows that the wave length 1 in free space , as defined by 
X = 2Tr/k, is 62.28 m. With speed of sound in warm air of 350 m/sec, 
the frequency f is f = kc/27r = 5-5T cycles /sec and the corre- 
sponding rotational speed of a compressor generating those waves is 
350 rpm which, again, is representative of industrial application. 

The argument of the Bessel functions will be, at the limits, 
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(kE^) = 0.02 and (kEg) = 0.037^+8 

This argvtment is << 1 and the wave length is much larger than the 
dimension of the bend which satisfies the requirements of the theo- 
retical treatment. 

4.'2 Bessel Functions and Characteristic Functions of the Problem 
4.21 Solution pertaining- to the real root . - Using the selected 
and calculated values of .‘.a and 'ws have, by equation (1*20) » 

the real root equal to 

V = 0.02828 (4.1) 

o 

tThe' evaluation of the Bessel functions of order . ±0 . 02828 for 

argimients (0.02), (0.03) and. (O. 037^89) and of the derivative of 

Bessel functions of order ±0.02828 for arguments kE^ = (0.02) and 

kEp = (0.037^8) is based on the expansion of J (kE) in terns of 

o 

increasing powers of the argument . 


1 

(±v)! 


(kr/2)^ 

:±v + l")T 



and on the expansion of (kr) given by equation (1.12). 

Because (kr) << 1 we may neglect the second and the following 
terms of the expansions. However, in order to evaluate the degree 
of approximation thus obtained we will keep the second teim of -the 
expansion. Since the expansions contain alternating terms, the error 
resulting from limiting expansion to a given term will be smaller 
than the first term which has been dropped. 


Without entering into details of algebra the folloving numbers 


are obtained, 

(0.02) = 0. 87789 [1.015T93 -- 0.0000987] = 0.891667 

o 

(0.03) = 0.8880llt[l. 015793 - 0.0002222] = O.9OI8U12 

o 

(0.037^8) = 0.89361 [1.015793 - 0.0003469] = 0.9074127 

o 

tr_ (0.02) = 1.13909 [0.983158 - 0.0001012] = 1.119790 

o 

(0.03) = 1.126108 [0.98315 8 - 0.0002276] =! 1.106886 

o 

J_ (0,03748) = 1.119056 [0.983158 -0,0003553] = 1.099811 

These numbers are shown on figure 10 . 

Keeping only the first two terms of the expansions (1.12) we 
obtain 

(0.02) = 1.25215 

o 

(0.03748) = 0.668134 
o 

(0.02)- = -1.59491 

o 

^°**^3748) = -0.851066 

~ o 

To check the results we substitute the values of the derivatives 
into equation (1.9) 

(kR^)j^^ (kRg) - (kR2)J_;_^ (kR^) = 0 

00 00 

The result is 

-1.06566 + 1.06561 = -0.00005 



(a) Bessel functions of order ±0.02828 in the 
interval of the argument (kr) from 0 02 to 0.03748. 



Figure 10. - Bessel functions of fractional and imaginary order. 


90 


Equation (1.9') is verified with an error of x 100 = 0.001 

percent. For the three radii the value of F (r) is 


= 0.2 


r = 0.3 


Eg = 0. 37^+8 


F = -2.82420 

V 

o 

F = -2.82434 
o 

F = -2,. 82437 

o 


with maximmn variation of F^ over to Eg of 2 ~ 8 ^ 2 8 ^ ' ” 

o 

0.006 percent. This result indicates that the wave is very nearly 


plane. Using the equation (2.4) which -was established without the 


second and following terms of expansion of J_^^(kr) we obtain 
F = -2.8240 

V / 

o 


The discrepancy between this value of F^ and the values eal- 

o 

culated using the first and the second terms of this expansion is 
negligible, as assumed in the evaluation of the roots of the wave 
equation. 

4.22 Solutions •pertaining to the imaginary roots. of -the problems . - 

To calculate the F^ pertaining to the imaginary roots , for 

m , 

m=l, 2, 3. . .»we start by evaluating the degree of approxi- 
mation of our mathematical theory. The pure imaginary roots are 
given by equation (1.17\) and only the first term of this equation 
was , considered in analysis because the second term was assumed small. 
Calculations confirm the theoretical development: 
m = 1 = i5 - iO. 000083 

m = 2 Vg = ilO - iO. 000040 

The second term is very small and decreases with increasing m. 
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To calculate the Bessel functions for v,m=sl, 2,3. . 
ve calculate functions G^(kr), F^(kr) and r(l±iy) for different 
y and (kr). 

We do. not intend to go into details of the algebra. However, 
it may be worthwhile to show the technique of work and to give a 
sample calculation. We select 
(kR^) = J^^(0.02) 

and propose to solve the equation (2.8a) with 
cos[5 In 0.02] = cos[-19.5601] = 0.75T997 
sin[5 In 0.02] = sin[-19 .56oi] = -0.652258 
Substituting 

0.757997 - iO. 652258 = 2 ^^ r(l + i5)J.^(0.02) (4.2) 

The- left-hand side of the equation is a complex number whose modtilus 
is 0.75799^ + 0.652258^ = 1 and the argument is 

2 ir - arctg = 5,57258. Substituting into (4.2) and taking 

0. (5 i99 ( 

logarithm we have 15.57258 = 13.^655 - 6.13032 + 13.815898 + 

In 'J.^(0.02) 
i5 

where 

-6.13032 + 13.81898 = In r(l +15) as given by tables of functions, 
reference (32). 

Finally 

6.13032 - il. 708818 = In J.,(0.02) 

i5 

or 

459.6(cos 1.708818 - i sin I.708818) = J .(0.02> and 

^ P’ 

J^^(0.02) = -63.236 - 1455.229 
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Calculating the Bessel functions for arguments (0.03) and 
= 0.037^89 we obtain 

J^’^(0.03) = 437-711j + il40.l3 
0.037489) = 63.236 + i455.229 

The Bessel functions (as 'approximated by limited expansions )' fo’a? 
the same three arguments but of negative order are complex con;3ug- 
ates of the above results. Calculating the derivatives of Bessel 
functions with respect to the argument, for arguments (kR^) 

0.02 and (kR^) = 0.037489, we obtain 
J^^(0.02 = 1.13807x10^ - 10.15805x10^ 

J|^( 0.037489) = -6.0712x10^ + 10.843174x10^ 
and the same functions but of negative order are the complex con- 
jugates of the above results . Again these conclusions are re- 
stricted to approximated expressions for the Bessel functions. 

The calcTilated functions are shown on figure 10. On this fig- 
ure we also indicated two particular- values of J^^(kr) namely, 
the pure imaginary J^^( 0.02043) = -i459.6 and the real . 
1^^(0.02815) = 459.6. -To check the calculated derivatives, we 
substitute the nmerical values of the functions into equation (1..9/). 
The results is 

J^^(0.02)J_|_^^(0. 037489) - 1^^(0.037489)1^^(0.02) = - i0.0231 - 0 

The accumulated error equals 0.12 percent and may be tolerated. 

4.3 ■ The Two Physical Systems 

The degree of accuracy of the derived basic 'functions has been 
established in the preceding section. The numerical application of 
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the integrated eq,uations for the infinite coil and for the 90 degree 
hend followed hy a straight infinite duct 'may now he done without 
going into details of the component functions . 

4.31 Infinite coil . - !Ehe tangential and radial vibrational velo- 
cities are given by equations (2.22) and (2.23). The-numerical appli- 
cation is straightforward. Ttie alternating series in the two equa- 
tions are converging rapidly as long as 6 is not zero. 

At _ 0 - ir/64, tt/8 and tt/ 4 the seventh, the- fifth, and the 

second term of the series, respectively, are already negligible. 

At 0 = 3ir/8 the first terms of the series was smaller than 0.1 
percent of the value of the simple wave term in equation (2,22) . 

At 0 = 0 the convergence is slow. Typically for 0 = 0 and at 
r = v(E^,0,t) = v^ 

and 

CO 

a - 1 p a cos(ni7r) - 1 „ 

In a In a(l + v^)' 

m=l 

and for 6=0 and r = 

f , 

00 

1 a. •- 1 ^ \ a - cos(m'ir ) _ ^ 

a In a a / ’ In a(l + y^) 

• ‘—r m 

m=l_ , 

With 30 terms taken into consideration the error was still of the 
order of 1 percent. The calciilated values of v(r,0,t) are plotted 
on figure ll. 

To better show the rapid velocity changes -near the inlet section 
of the bend and. the important final velocity dependence ,(l/r) on the 




(a) Tangential vibrational velocity distribution near inlet of an infinite coil. 




(b) Magnitude and direction of vibrational velocities at section 0 = 7r/I6, 
Figure 11. - Infinite coil. Distribution of vibrational velocities. 
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radius, the calculated values are superposed. On the abscissa is 
the velocity magnitude in terms of v^. 

About 25 percent of the velocity redistribution take place be- 
tween 9 = 0 and tt/64. By 9 = tt/8 the distribution character- 
istic to potential vortex is nearly completed and changes become 
slow and it takes, from then on, an angle Tr/4 to reach the steady- 
state distribution. 

The decaying radial vibrational velocities were calciilated for 
one cross section, at 0 = it/i6, and are compared with the permanent 
radial vibrations. The values are 

Radius Decaying ' Permanent 

u(r,ir/l6,t ) u(f,iT/l6^t)' 

0 ; 0 

R, + (R^.- Rj/h -0.0860 V -0.002475 V 

R, + (R^ - R, )/2 -0.0884 V -0.002636 V 

± id 1 >0 .0 

R, +-3(Ro - R-,)/4 -0.0480 V , -0.001635 V 

1 d 1 ' . O '<0 

Eg *0 0 

The permanent radial oscillations are one. order of magnitude smaller 
than the decaying oscillations at ir/l6. By 'Vectorial addition, the 
magnitude and direction of vibrational velocities may be obtained. 
For section 6 = it/16, figure 11 shows that the wave profile is not 
straight. It approximates in the vicinity of this section, a cosius 
curve. Beyond section 0 = Sv/S the radxal velocities are so small 
that the wave profile is very nearly radial. 
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4.32 90 degree 1)6113. foilowed by a straight duct >. - The wave pattern 
in a tend followed hy a straight .duct is much more involved. than the 
wave propagation in an infinite coil. So is also the numerical 
application. 

The problem of evaluation of integrals A^. was already dis- 
cussed. and is reported in appendix 4. The evaluation of coeffic- 
ients D or rather of the ratio (D /iC. ) j given by matrix (2. ,32a), 
m mm 

was simplified, by elimination of small terms and reduced to matrix 
equation (2.32b) in This matrix, in principle infinite was , in 

the present numerical calculations, arbitrarily -limited to six rows 
and columns. It is as follows: 


_ , 0^ 0.3288 0.401^, Q.105J g . 0.978 , 

2 - 1.9935 Bg + 10^ ^3 lo"^ ^4 ^ 10^ ^5 ^6 

_2 = — 0 208 6 ' - 2 001 6 ' — 6 * + ^ ‘ g t _ 0 « 1358 _ , 0 . 2 T 4'3 - g , ■ 

-2 - - 0.200 - 2.001 - “io^ ^3 + 10 ^ ^4 10 ^^ ^5 ^ 10 ^ ^6 


_2 = 


-2 = 


0 7L5 «• 0.^277 p O non R' 0.4887 g. . 0.l4l 0.786 

- Bg - 2.000 B3 - -3;^- Bi^ + B. - B 


- 0.4128 B{ + B^ - B^ - 1^995 B^ - 3^ 


2- 0.455 g. ' 0.120 0.181 . 0.108 pg 0.314 . 

" 10 ^1- 10 ^2'^ ^0^ ^3 10 ^4 ^ ^5 “ ^10^ ^6 

-2 = - 0.7615 Bj^ + - - ^ 3^ B^ - 0.2699 B^ + 0.135 B]^ - 0.844 B^ - 1.5768 B^ 


By inspection, we see that the first four B^^'s will be close to 
iinity. To have precise values, the matrix was solved by Cramer's 
rule on a digital computer. ’Furthermore, in order to ascertain that 
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the arbitrary limitations of an infinite matrix to a square matrix 
6 x 6 will not result in unrealistic values of five computer 

runs were requested to obtain solutions of matrices 2 x 2 , 3 x 3 and 
6 x 6 . The result is 


Approximation 


^2 


^3 






^6 


1st 1.00326 

2 nd 0.99998 
3 rd 0.99896 
4 th 0.99894 
5 th . 0.99894 
6 th . 0.99893 


0.89374 



0.89374 

0,62598 


0.89379 

0.62579 

0 . 

0.89378 

0.62650 

0 . 

0.89381 

0.62640 

0 . 


T 9069 

7 T 137 1.00425 
77365 1.00024 0.24667 


The variation of coefficients 6’ over the- range of terms taken 

m 

into consideration is negligible. Consequently, we assume that ar- 
bitrary limitation. of the matrix to 6x6 is permissible. The values 
obtained by -the sixth .approximation will be used in calculations. 

The calculation of -the tangential and radial vibrational velo- 
cities [eqs. (2.35) and (2. 36 )] is direct. Figure 12 gives the re- 
sults of the .numerical application. At 6 = 15ir /32 -the distribu- 
tion of tangential velocities is no more inversely proportional to 
radius. The wave is. deformed, by significant, positive radial 
velocities. The series in the equations of motion are still limited, 
however, to one significant term only. At • 9 » 3111/64, there are 
two significant terms in the series. We conclude that in^ortant 
changes in the values of vibrational velocities are confined to the 
last tt/32 of the it/2 bend. 

By the bend outlet the wave profile is no more radial. The 
effective vibrational velocity equal + u^ 


is in r positive 




12 1.1 1 .9 .8 


Vo 

(c) Straight duct inlet Distribution of axial vibrational vel«Ities. (di Vibrational velocities at o. 25 

2L 

Figure 12. - 90“ Bend followed by a straight duct Distribution of vibrational velocities. 
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direction, except at 'bo'undaries where it is tangential. This is al- 
so shown on figure 12. The direction of the effective vibrational 
velocities gives the wave front orientation. 

The calculation of the vibrational velocities in-the straight 
duct that follows the 90 degree bend is straightforward. Results of 
application of equations (2 ►39') and (2. ■■40) are also reported on fig- 
ure 12 . 

In the straight duct the waves are strai^tening out relatively 
fast. By x/2L = 1 there is already an almost straight wave simi- 
lar to that which was generated at the bend inlet. Essentially, the 
process of straightening out of the wave is confined to the straight 
duct. Figure 13 shows that the tangential velocity "v = remains 
in the bend in its equilibrium position at approximately - “ 

l)/ In a, which in this application, is at r/R^ = 1 . 392 , and only 
in the straight duct gradually bends closer to the center line of 
the duct. By x/2L =1.5 the tangential vibrational velocity (v^ - 
6 ^ - 6^) is at the center line and, everywhere. 

The verification of the simultaneous solution of the equations 
of motion for the bend and for the straight duct is shown on figure 
lU. The calculated tangential vibrational velocities for r = R^ 
and r = R^ for positions starting at the vibrating piston, through 
the bend and in the straight duct are taken from figures 11 and 12. 
The rapid changes in velocities by the bend exit are well illustra- 
ted. In spite of this, the values .calculated for the bend match 
very well the values calculated for the straight duct. At x/2L = 1 




0 . 25 . 50 . 75 1.00 1.25 1.50 

Nondi men siona I distance, x/2L 

Figure 13. - Propagation in bend-straight duct system. Changes 
in location of the tangential vibrational velocity Vq. 



0 .25 .50 .75 1.00 

X/2L 


Figure 14. - Propagation in bend-straight duct system. Tangential vibrational 
velocities at bend’s curved walls, at r = and r = R 2 . 
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the vihrational velocities at the walls of the straight duct are, 
within a fraction of a percent , eq.ual to v^. 

For comparison purposes, the two tangential lines shown indi- 
cate the values of the vihrational velocities for an infinite hend. 
This graph should be compared with the well known and often pub- 
lished curves for potential fluid flow in elbows . Those curves give 
the fluid velocities of the streamlines at elbow walls , as estimated 
by semi- graphical solutions. The calculated curves on figure l4 
are similar. -The derived deviations from a logarithmic distribution 
of velocities in a bend are in the case of long wave mot-ioir, -too 
small to be noticeable. 



'5.0 CON CLIjDING EEMAEKS MD RECOMMENDATIONS 


Propagation of long vaves in a two-dimensional system has heen 
analyzed. The -acoustic approximation has heen used and the Helmholz 
equation was integrated for two sets of boundary conditions. The 
two physical systems taken into consideration are : an infinite bend 

approximating a coil and a 90 “^ 'bend followed by .a strai^t , infinite 
duct approximating a -typical indistrial piping system. 

Implicitly the analysis gives also a solution for a system con- 
sisting of an infinite duct connected to a 90 ° elbow followed by a 
second infinite duct. The results of the analysis indicate that 
bending of a straight duct profoundly modifies the propagation of 
waves in that duct. The bend engenders the following phenomena: 

a) a set of attenuated axial and radial waves. which modify 
the plane wave generated at the duct's inlet (figs. 5 and 
6 ). 

b) a continuous radial, standing wave whose radial vibrations 
are . sustained by the curvature of the bend. In a straight 
duet these vibrations would be quickly attenuated (fig. ?)• 

c) a vortex-type distribution of the tangential vibrational 
velocities (fig. 3). 
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d) an increase in the phase velocity which is proportional to 
the sharpness of the bend (fig. 4), 

e) a decrease of the amplitude of the transmitted wave. The 
radial reflections at the bend make the transmitted wave 
weaker than it was before. 

The analytical results have been qualitatively verified by data ob- 

(25 ) 

tamed by Lippert^ . In spite of the differeiices between the sys- 
tem considered in the analysis and the system which has been used in 
experiments there is an agreement on the trend of the change of 
parameters of the waves crossing those bends. 

The presented analysis is not directly applicable to real flows 
because it is based on a linearized equation valid for acoustical 
waves in stationary medium only. It gives, however, an idealized 
picture of wave's behavior in bends, in general. This work extab- 
lishes a method to estimate capability of bends to transmit waves. 

It may be useful in the design of rigid accoustical duct systems. 

It also opens possibility to calculate transmittability of non-rigid 
bends which induce large but as yet not calculable phase changes in 
the progressing waves. 

The obtained results should also be useful in advancing the 
theory of waveguides in electromagnetic application. In spite of 
the fact that the electromagnetic theory is concerned with high fre- 
quencies only, the developed reactions, particularly for the decay- 
ing field of velocities, should be of value. In the field of low 
frequency oscillations, analysis of the observed ”Pogo" 
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instabilitzes of rocket systems should henefit from the developed 
eq.uatlons , 'While older rochet designs were provided with straight 
feed piping systems, the newer designs will have bends and elbows 
which must be analyzed. 

A more complete research may now be suggested to get a better 
understanding of the phenomena associated with the propagation of 
pulsating flows in curved pipes. From the stability of flow view- 
point the most interesting case to investigate would be a free pip- 
ing system, (i.e.j a bend in a piping system with one or two degrees 
of freedom). Across the elbow an important phase change should be 
obtained. From the heat transfer point of view it might be inter- 
esting to solve a tri-dimensional Helmholz eq.uation and wave propa- 
gation combined with fluid flow. 
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APPENDIX 1 


The Equation of Continuity 

¥e define dilation 6 =(u “ ^o^^o ^ is an elementaiy 

volume and the subscript zero refers to the initial state'. If ^ 

1 

are displacements in three-dimensional space 


3x Sy - 3z 


The coefficient of compressibility of the fluid under consideration 
is defined by' 

u - u 
o 

X = - 

hence 

P = P - - and - - div V (Al.l) 

o X X 

■where V is the velocity vector. Flow of material out of the ele- 
mentary volume reduces the pressure of the compressible fluid. The 
rate of change of excess pressure is proportional to the divergence 
of V. For incompressible fluid div V = 0. 

For an isentropic process pu^ = const and -udp + yp^U = 0 
and considering the definition of x 

dU _ ^ 

^ ~ ~ udp " pY 



If we put Y ^ where c is. the velocity of sound, as 


established by Laplace in l8l6, we have 



(A1.2) 
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APPENDIX 2 


The Orthogonality Condition 


Fxmotion (kr) ^ v 

m m m mm 

every must satisfy the Bessel Equation. Written for two differ- 

ent Vj and *the two relations are; 


1 / 

p” +— F' + (k^ - -4 I P = 0 


F" + - F' + (k^ - -| F , =0 

^2 ^ ^2 \ ^ / ^2 


Multiplied hy , rF ^ »and rF . respectively,. integrating, 
and rememhering that 


F^ .(^^2.^ “ ^^^2^ = 0, we get 

m m 



and a similar equation for the second expression.. Subtracting, we 
obtain 



or 
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and. assuming that 


Vg we detemllne -the orthogonality condition 



(A2.1) 


From the above analysis we also conclude that no complex roots 
can exist. For if were such a complex root, its conjugate v 
would also be 'a root and then F_ = F and the product F F_ would 
be necessarily a positive number, different from zero.. Consequently. 

must have either a ptare imaginary value or real va:|.ues for other- 
wise 


Vg - = (v + v)(v - v) = 4EevImv 

would not vanish. 

It should be pointed out that a similar relation but for Bessel 
functions of pure imaginary argument (ikr) has been discussed by 

(33) 

Bocher'"^-^-' . 
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APPENDIX 3 

The Integral ar 

Function F = J (kr)J' (kR, ) - J (kr)J’(kR-.) satisfies the 

V V -V 1 -V V 1 

Bessel eq.uation 


F" _ i. F' + (k^ - F = 0 
V r V r / V 


(A3.1) 


Let us differentiate (A3.l) ‘with rgspect to v and put 
9F 


(v) 


V 

3v 


FV V + - FI V + k^I 
(v) r . (v) 


(v) 


-^F, v-^F = 0 
(v) V 


(A3. 2) 


Multiplying (A3.l) By (A3. 2) By 'we have 

J ^Vfv) * f ^(vfv 


1 

r "(vrv‘ 




- 2v / - F^ dr = 0 

■ V 


rF^y^Fj dr +| F^F^^^dr + k^y rF^F^^^dr 




F F/ ,, dr = 0 
V (v) 


SuBtractlng and integrating By parts the integrals containing 
the second derivatives of F, we have 


rF FI ,, 
V (v) 




R„ / 

2 


2 / 


rF, vF* 
(v) V 

T, - 2v/ 

?1 

J 


Applying the Boundary condition F^(kR^) = F^(kR2) =-0 we have . 
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-f 2 dr 
r V • 


— rF f; n 
2v V (v) R, 


8J'(kr) 

V 




9J' (HE ) 3J' (Er) 3J'(kS ) 

* - ^17— — TT-' 
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APPEMDIX k 


The Integral A = 
inn 


cos 







r 

cos 

^ 5“ 
m 


dr 


The symbol stands for the integral 



'R. 


cos 


'R. 


■) 


ijT - E^) 


cos 


r 

y In — 
m R 


U 


dr 


which is part of integral 




Rmn 


R, 


fn< 


F (r) cosk (r - R ) 
n 1 

m. 


dr, m ^ 0. 


The integral and, ipso facto, the integral appear in 

eq.uations (l.&O) and (l.4l). The integral A^ has no simple solu- 
tion; however, the integral A^^ will he easily obtained when the 
basic integration of- A^ will he performed. 

The problem is thus reduced to evaluation of A . The pro- 
• mn 

cess of integration of the definite integral A^ begins with trig- 
onometric transformation to split the integrand into four definite 

integrals of products of sine and cose of (? r) and p In r. Next, 

n m 

expanding cos C r and sin 5 r in series we obtain four alter- 
n n 

nating series containing integrable terms of the types 



cos 

cos 


y In'-r 
m 

y In r 
m 



11 ^ 


Using the definite integrals developed earlier, eq.nation ( ’), 

and applying the integration limits, after simplification, ve 
finally obtain 


A = B- cos f? R^V-Tp-T — ^ [a cos(nnr) - l] 
mn "1 VnlMl + u^ 


n 

3R2 

21 

(3^ >‘5) 

n 

5Sj_ 

iT 

P" * 

\ 



Ur| 

31 


n 


51 



- 1 ] 


•where m=l,2,3... n = l,2,3.. 


U = m 


C ~ n 


m In a ’n 


a = 


h 


For the numerical application ir/ln a = 5, 

Eg - Ej_ = 0.1T489m, E^ = 0.2m 

The derived series solution for A is difficult to handle . 

mn 

The two converging series exhibit behavior similar to that of the 
series of sine and cosine of argument »1, The second and several 



115 


following terms of the series grow very fast and attain very high 
values hefore the factorials in the denominators force the decrease 
of the terms. 

An attempt was, made to compute on a digital computer, 

IBM 709^ Mod II. Good results were obtained for values of m = 1, 

2, 3 . . ♦ 6 and for n = 1 and. 2. For n = 3 the results were 
erratic and for n > 3 even the double precision was not suffic- 
ient to carry out calculation and completely erroneous results were 
printed out. 

A second approach consisted in use of an existing routine, 
recorded on tape,. of -integration by Simpson’s -Rule. A con^ilete 
set of values was,, obtained for n and m up to 6 and in order to 
study the behavior of the function several values for higher’ and 
m. All calculated are reported on figure 15. The complete 

set of results is presented in the form of a matrix: 


Values of A 

mn 



n = 1 

n = 2 

n = 3 

n-= 4 

n-= 5 

n = 6 

m — 1 

0.08583 

0 . 0104l 

0.00140 

0.000509 

0.000148 

0.0000849 

m - 2 

-0,02263 

0.08187 

0,02076 

0.004048 

0.00144 

0.000466 

m - 3 

0.00324 

-0.03208 

0.07523 

0.02976 

0.00762 

0.00276 

m = 4 

-0.00232 

0.00667 

-0 .o 4 o 69 

0.06633 

0.03715 

0.01189 

m - 5 

0.00054 

-0.00319 

0.01109 

-0.04774 

0.05559 

0,04262 

m “ 6 

-0.00079T 

0.000934 

-0.00463 

0.01620 

-0.05285 

0.04350 


Along the principal diagonal A^ steadily decrease. Along 

the' super and sub-diagonals the values of A increase, with A 

" mn mn 
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for m = n varying- along , two sine curves. By m = -n = 6 the 
matrix is approximately tridiagonal. Figure 15 shows ■ a. series of 
functions A^f(m) which are rapidly attenuated for n < 3 but ex- 
hibit large oscillations for higher values of n. The values cal- 
ciilated by Simpson’s Eule agree well with values obtained by the 
series solution for n = 1 and 2. 

f 

For n 5= 1 and m = 1 . . . 6 we have 

m 1 2 3 4 5 6 

Series 0.08583 -0.02263 0.00324 -0.00232 0.00054 -0.000797 

solution 

Simpson's 0.08577 -0.02266 0.00327 -0.00234 0.00055 -0.000804 

Eule 







